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INTRODUCTION: 

We  recognize  that  many  patients  with  Multiple  Sclerosis  (MS)  suffer  from  cognitive  impairment  at  some  point 
in  their  disease  course.  However,  characterization  cognitive  change  in  patients  with  MS  has  been  difficult  to 
pinpoint,  and  is  hampered  by  poor  quantitative  markers.  We  have  two  hypotheses:  1)  conventional  imaging  is 
insensitive  to  gray  matter  (GM)  changes  known  to  exist  in  patients  with  MS,  and  2)  ultra-high  MRI  field 
strengths  (7T)  would  allow  an  opportunity  to  study  the  myelination  and  metabolic  changes  of  the  cortical  GM  in 
patients  with  MS  and  known  cognitive  impairment.  The  purpose  of  this  proposal  is  to  develop  and  implement  a 
targeted  quantitative  magnetization  transfer  (qMT)  and  chemical  exchange  saturation  transfer  (CEST)  MRI 
imaging  paradigm  at  7T  to  detect  and  quantify  the  level  of  myelin  loss  (qMT),  protein/peptide  changes  (amide 
proton  transfer  CEST),  neurotransmitter  deficiencies  (GluCEST)  in  the  GM  of  patients  with  MS,  and  to  relate 
these  findings  to  neuropsychiatric  evaluation  outside  the  MRI  scanner.  The  scope  is  to:  1)  develop  novel,  high- 
resolution,  high  field,  quantitative  MRI  methods  sensitive  to  myelination  and  neurochemicals  for 
implementation  in  the  cortical  GM  of  human  populations,  2)  deploy  these  methods  in  patients  with  MS,  3) 
relate  these  findings  to  measures  of  cognitive  impairment,  and  4)  develop  a  lower  MRI  field  strength  alternative 
for  direct  patient  impact. 

KEYWORDS: 

•  Magnetic  Resonance  Imaging  (MRI) 

•  7  Tesla  (7T) 

•  Chemical  Exchange  Saturation  Transfer  (CEST) 

•  Magnetization  Transfer  (MT) 

•  Brain 

•  Cortical  Gray  Matter  (cGM) 

•  Multiple  Sclerosis  (MS) 

•  Functional  MRI  (fMRI) 

•  Pool  Size  Ratio  (PSR) 

•  Amide  Proton  Transfer  (APT) 

•  Glutamate  (Glu) 

•  Myoinositol  (ml) 

•  Cognitive  Impairment 

OVERALL  PROJECT  SUMMARY 

Task  1.  IRB  Preparation  and  Human  Subjects  Approvals. 

Completed 

Task  2.  Develop,  optimize  and  implement  advanced,  quantitative  Magnetization  Transfer  (MT)  and  Chemical 
Exchange  Saturation  Transfer  (CEST)  in  phantoms  and  evaluate  minimum  achievable  resolution  and  the 
associated  reliability  of  derived  indices 

Simulation/Phantom  Studies 

The  objective  of  this  task  was  to  develop  a  best-practice  MT  and  CEST  acquisition  scheme  to  be  deployed  in 
Task  3  in  healthy  controls. 

Summary  of  Results/Progress  and  Accomplishments 
CEST- 

APT-CEST  -  We  have  optimized  through  simulation  and  phantom  studies  a  single-power,  whole  brain  APT 
CEST  acquisition  for  deployment  in  healthy  volunteers  and  patients  (Y2Q1).  We  began  with  the  protocol 
presented  by  Jones  et  al  (1)  for  whole  brain  coverage  at  7T  and  increased  the  in-plane  resolution,  the  coverage, 
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and  modified  the  saturation  scheme,  readout,  and  fat-saturation  (now  a  binomial  excitation  pulse)  to  minimize 
distortions  in  the  phantoms  while  maintaining  sufficient  contrast  to  noise  (CNR)  for  the  APT  CEST  signatures 
in  a  reasonable  scan  time  (9  minutes).  The  relevant  scan  parameters  are  as  follows: 

•  Whole  Brain  (33  slices)  3D  Gradient  Echo  (FFE)  with  multi-shot  EPI  (factor  =  7)  readout 

•  1.5  x  1.5  x  5mm3  acquired  resolution 

•  CEST  RF  Saturation  B 1  =  2[jT  x  25ms  (each) 

•  64  offset  frequencies  (Dw  =  -5  -  5ppm,  Aw  step  =  0.2ppm  +  14  no  Saturation  Scans) 

•  Total  Scan  Time  =  9:10 

The  simulations  were  performed  using  a  3  pool  model  of  the  CEST  effect  presented  by  Zaiss  et  al  (2)  and 
inputting  estimates  for  Tl,  T2,  pool  sizes  and  exchange  rates  for  the  macromolecular,  labile  (amide  protons) 
protons,  and  water  as  follows:  Simulations  were  designed  to  identify  optimal  CEST  preparation  (RF  irradiation 
power  and  bandwidth)  parameters  for  APT  contrast.  All  simulations  were  carried  out  utilizing  the  scripting 
environment  in  MATLAB  2012b  (Mathworks,  Natick,  Massachusetts)  on  an  Apple  iMac  (Cupertino,  CA;  3.0 
GHz,  dual  core  CPU).  Theoretical  saturation  was  modeled  according  to  the  Bloch  equations  for  three  pools: 
bulk  water  (free),  semisolid  macromolecular  (conventional  MT),  and  mobile  macromolecular  (CEST)  pools. 
This  was  achieved  using  the  simple  matrix  solution  to  numerically  solve  the  Bloch  equations  (3).  This  model 
assumes  a  T2-dependent  Super-Lorentzian  absorption  lineshape  for  the  macromolecules  (4).  Physical  values  of 
exchange  rate,  Tl,  and  T2,  were  fixed  according  to  (5,6).  The  bulk  water  was  modeled  as  a  Lorentizan  with 
T1/T2  =  1538  ms/45ms.  The  semisolid  macromolecular  pool  was  modeled  with  Tl/T2/exchange  rate  = 
1600ms/0.01ms/20  Hz  with  offset  =  -2.34ppm  (7).  The  labile  proton  pool  was  modeled  as  a  Lorentzian  with 
Tl/T2/exchange  rate  =  500ms/20ms/50Hz  with  offset  =  3.5ppm  and  a  concentration  of  0.001  compared  to  bulk 
water  (1.0)  and  macromolecular  pool  (0.1).  Therefore  our  simulations  included  a  single  RF  irradiation,  a  brief 
delay  for  spoiling,  an  on-resonance  excitation,  and  a  delay  for  readout.  The  Bi  amplitude  (power)  was  varied 
over  amplitudes  from  1  pT  to  3  pT  while  holding  the  duration  constant  at  25  ms.  The  Bi  amplitude  was 
subsequently  fixed  to  1  pT  while  the  pulse  duration  was  varied  from  0  ms  to  60  ms. 

For  phantom  scans,  we  performed  the  above  paradigm  at  various  in-plane  resolutions  as  low  as  0.75mm2, 
but  determined  through  curve  analysis  that  at  these  resolutions  the  CNR  was  insufficient  to  parse  out  the  CEST 
effect  from  the  background  noise.  From  these  simulations  and  phantom  studies,  the  above  pulse  sequence 
paradigm  was  chosen  to  maximize  CEST  contrast  derived  from  APT. 

GluCEST  -  For  glutamate,  we  performed  the  same  simulations  as  above,  but  we  modeled  the  off-resonance 
saturation  as  given  in  (8),  and  assumed  an  exchange  rate  for  glutamate  =100  Hz,  pool  size  =  0.001  compared  to 
bulk  water  (1.0)  and  macromolecular  concentration  (0.01)  and  Aoj  =  3.0ppm.  For  the  GluCEST  acquisition,  due 
the  concern  of  overlapping  resonances  (GABA,  Glutamate,  and  other  Amines)  a  high-spectral  resolution 
acquisition  needed  to  be  obtained  (Au)  =  0.2ppm  spacing),  thus  scan  time  becomes  prohibitive  for  extremely 
high  resolution.  However,  it  should  be  pointed  out  that  the  amine  resonances  that  may  reside  in  juxtacortical 
WM  will  be  significantly  less  than  the  adjacent  GM,  so  a  slightly  poorer  resolution  acquisition  will  not  be 
problematic  if  a  high  spectral  resolution  scan  is  obtained.  We  therefore,  decided  to  utilize  a  scan  very  similar  to 
that  which  as  been  presented  by  Dr.  Reddy  (8,9)  and  thus  we  will  implement  a  single  slice  GluCEST  acquisition 
in  vivo  at  a  resolution  of  1.9  x  1.9  x  5mm3. 

•  Single  Slice  (2D)  Gradient  Echo  (FFE)  with  multi-shot  TFE  (40  shots)  readout 

•  1 .9  x  1 .9  x  5mm3  acquired  resolution 

•  CEST  RF  Saturation  BI  =  4.25  pT  x  10ms  (each)  x  lOOsegments  at  90%  duty  cycle 

•  50  offset  frequencies  (Aw  =  -5  -  5ppm,  Aw  step  =  0.2ppm) 

•  Total  Scan  Time  =  1 1 :36 


qMT  - 

We  have  chosen  to  implement  the  selective  inversion  recovery  (SIR)  quantitative  MT  (qMT)  to  quantitatively 
extract  the  pool  size  ratio  (PSR),  which  has  been  shown  to  be  reflective  of  myelin.  We  have  developed  the  SIR 
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approach  at  7T  as  discussed  in  the  Q1Y1  and  Q2Y1  progress  reports.  This  pulse  sequence  has  been  shown  in 
previous  reports,  published  (10)  and  provided  in  Appendix  3.  However,  in  Q3Y1,  we  studied  via  simulation  the 
impact  of  partial  volume  effects  where  we  know  there  is  a  non-negligible  MT  effect  in  GM  and  certainly  a 
strong  effect  in  WM.  Thus,  in  juxtacortical  voxels  where  a  blend  of  GM  and  WM  may  occur,  poorer  in-plane 
resolution  results  in  an  inability  for  the  model  to  remain  stable  when  deriving  the  PSR  values  (i.e.  two  different 
PSR  values  may  fit  equally  well  when  there  are  two  populations  within  a  voxel).  Therefore,  we  proposed  a 
reduced  number  of  slices  but  increased  the  in-plane  resolution.  From  simulations,  we  feel  that  this  provides  the 
most  robust  acquisition  method  to  be  deployed  in  patients.  Thus,  from  our  phantom  studies,  we  have 
determined  that  a  1  x  1  x  2mm3  acquisition  with  5  slices  sampled  at  14  TIs  (TI  =  6ms,  10ms,  16ms,  26ms, 
42ms,  68ms,  110ms,  178ms,  288ms,  468ms,  760ms,  1233ms,  2000ms,  8000ms)  will  be  performed  in  patients 
with  MS  and  healthy  controls. 

•  Inversion  prepared  3D  Gradient  Echo  (FFE)  with  multi-shot  TFE  (2  shots)  readout 

•  1  x  1  x  2mm3  acquired  resolution  (5  slices) 

•  14  Inversion  Times  (TI  =  6ms,  10ms,  16ms,  26ms,  42ms,  68ms,  1 10ms,  178ms,  288ms,  468ms,  760ms, 
1233ms,  2000ms,  8000ms)  at  a  constant  delay  time  (TD  =  2500ms) 

•  Total  Scan  Time  =  10:11 


Conclusion  of  Task  2:  Simulation  and  Phantom-optimized  qMT  and  CEST  acquisitions 

Through  simulation  and  phantom  studies,  we  have  devised  a  final  protocol  to  be  deployed  in  healthy  volunteers 
and  patients  with  MS.  A  summary  of  the  protocol  is  given  below,  and  a  complete  protocol  is  given  in  Appendix 
1. 

Final  Summary  of  Protocol  implemented  in  healthy  controls  (Task  3)  and  patients  with  MS  (Task  5). 

•  Constant  RF  APT  CEST  -9:10 

•  Constant  RF  GluCEST  -11:36 

•  SIR  qMT  -10:11 

•  Bloch-Siegert  B 1  mapping  -  1 :42 

•  Dual-echo  B0  mapping  -  :04 

•  T 1  w  MPRAGE  Anatomical  -2:12 

•  fMRI  Resting  State  -  8:34 

•  fMRI  N-Back  task-  8:30 

•  fMRI  Trailmaking  task  -4:14 

The  current  scan  time  for  all  scans  is  approximately  1  hour. 


Task  3  —  Implement  current  best  practice  for  MT  and  CEST  in  healthy  volunteers  and  evaluate  reliability 

The  objective  of  this  task  was  to  implement  a  best-practice  MT  and  CEST  acquisition  scheme  in  healthy 
controls. 


Summary  of  Results/Progress  and  Accomplishments 

The  above  protocol  has  been  implemented  in  20  healthy  volunteers  at  the  close  of  year  1 .  We  have  additionally 
repeated  this  paradigm  in  8  healthy  volunteers.  We  have  further  scheduled  the  remainder  of  the  healthy  controls 
to  be  scanned  in  the  coming  month.  There  is  one  delay  to  report  in  that  at  the  close  of  Year  1,  our  SOW  stated 
that  we  would  have  recruited  50  healthy  volunteers  into  the  study.  As  pointed  out  in  quarterly  reports  Q1Y1 
and  Q2Y1,  we  struggled  initially  making  the  phantoms  to  study  the  impact  of  resolution  on  the  final  protocol. 
This  resulted  in  less  than  the  expected  50  healthy  controls.  However,  we  have  already  scheduled  these 
remainder  healthy  volunteers,  and  have  a  40%  return  rate  on  for  repeat  visits  to  understand  the  reproducibility. 
We  will  complete  the  healthy  volunteers  and  repeat  visits  in  Q1Y1 
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In  a  follow-up  to  the  Q3Y1  report,  we  have  added  three  fMRI  scans  in  collaboration  with  Dr.  Paul  Newhouse, 
our  neurocognition  expert.  It  is  important  that  we  note  that  this  does  not  change  the  scope,  but  rather  offers  a 
unique  opportunity  to  study  the  cognitive  function  in  the  MRI  in  healthy  volunteers  and  patients  for  greater 
understanding  of  the  relationship  between  the  advanced,  quantitative  measures  and  outside-scanner 
neurocognitive  battery.  This  is  exceptionally  unique  as  neither  of  these  three  fMRI  scans  have  been  studied  in 
MS  patients  with  known  cognitive  impairment  and  creates  an  exceptionally  rich  data  set  to  mine  for 
understanding  neurocognitive  decline  in  MS  patients. 

In  20  healthy  volunteers  (with  8  repeat  acquisitions),  we  have  obtained  the  entire  proposed  MRI  protocol  as 
given  above.  Preliminary  results  follow  under  Task  4. 

Task  4  -  Analyze  the  derived  indices  in  healthy  volunteers  and  evaluate  reproducibility  (1  month) 

The  objective  of  this  task  was  to  1)  develop  an  analysis  pipeline  for  constructing  maps  and  deriving  indices 
reflective  of  GM  and  juxtacortical  WM from  the  quantitative  MRI  acquisitions  prepared  in  Task  3,  and  2)  to 
ascertain  these  indices  in  preparation  for  analysis  of  reproducibility. 

Summary  of  Results/Progress  and  Accomplishments 

Develop  an  analysis  pipeline  for  routine  analysis  of  data  generated. 

As  the  first  part  of  Task  4,  we  well  understood  the  need  to  1)  correct  patient  motion  in  an  individual  scan 
(motion-correction),  2)  co-register  data  across  scans  into  the  same  space  for  robust  analysis  (co-registration), 
and  3)  segmentation  of  WM  and  GM  for  histogram  and  descriptive  statistics  of  each  derived  index. 


WM/GM  segmentation  was  performed  in  FAST  using  3  classes  as  implemented  in  the  FSL  toolbox  (FMRIB, 
Oxford,  UK).  The  co-registration  was  performed  using  FNIRT  (non-linear  registration,  FSL,  FMRIB,  Oxford, 

UK)  to  put  the  quantitative  maps  into  the  MPRAGE  space 
such  that  the  segmentation  can  be  applied.  To  that  end, 
Figure  1,  shows  two  slices  of  an  APT-CEST  acquisition 
motion  corrected  and  co-registered.  The  magenta  color 
indicates  the  WM,  and  the  green  indicates  the  GM  derived 
from  the  MPRAGE  anatomical  acquisition  and  overlaid  on 
the  APT  map.  It  can  clearly  be  seen  that  the  agreement 
between  the  MPRAGE  and  the  APT  maps  is  high  and  the 
WM  and  GM  clearly  seen.  The  bottom  panels  show  the 
process  of  joining  the  MPRAGE,  the  WM/GM 
Segmentation  and  the  GluCEST- weighted  acquisition  for 
completeness. 


Figure  1  -  (Top)  coregistration  and  segmentation 
results  applied  to  APT  maps,  and  (bottom)  the  target 
(MPRAGE),  segmented  map,  and  GluCEST- 
weighted  acquisition. 


APT  CEST  analysis  and  results 

We  have  constructed  the  APT  CEST  maps  in  the  following  manner.  First,  the  CEST  spectrum  for  each  voxel  is 
normalized,  corrected  for  B1  drift  and  fit  to  a  single-lorentizian  (11)  and  the  minimum  spectral  intensity  is 
shifted  to  an  offset  (Aw)  =  0  for  BO  correction.  After  this  correction,  the  difference  between  the  data  and  the  fit 
create  a  Lorentzian  residual.  The  residual  between  Aw  =  3.25  and  3.75  are  integrated  and  termed  the  APT 
Lorentzian.  An  example  of  this  is  shown  in  Figure  2A  left  panel.  To  assess  reproducibility  we  created  a  single 
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CEST  spectrum  for  the  whole  brain  and  compared  visit  1  to  visit  2.  Figure  2A  right  panel  shows  the  results  of 
two  healthy  repeated  studies,  where  the  green/blue  curves  are  subject  1,  and  magenta/black  curves  are  subject 
two  over  two  times  points.  As  it  can  be  seen,  the  reproducibility  is  high  over  the  whole  brain.  Individual 
structure  assessments  are  ongoing.  It  can  also  be  appreciated  the  spectral  quality  of  the  CEST  spectrum  using 
this  analysis  approach. 
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Once  the  data  were 
segmented,  we  compared  the 
GM  averaged  over  all  healthy 
volunteers  (n  =  20)  and  1  MS 
patient  with  clinically 
diagnosed  cognitive 

impairment  via  histogram 
analysis  shown  in  Figure  2B. 
We  performed  the  same 
calculation,  but  also  examined 
resonances  at  2ppm  (hydroxyl 
and  sensitive  to  myoinositol) 
and  3ppm  (amine  protons 
sensitive  to  GABA  and 
Glutamate)  and  3.5ppm 
(amide  proton  transfer  - 
sensitive  to  pH  and  protein 
concentration).  It  should  be 
pointed  out  that  the  1  MS 
patient  examined  here  is 
actually  part  of  Task  5,  but  it 
is  important  to  show  here  as  it 
points  towards  the  sensitivity 
of  the  measurement.  In  this 
one  patient,  there  is  an 
obvious  downward  shift  of  the 
GM  histograms  at  2ppm  and 
3.5ppm  giving  the  impression 
that  we  are  detecting  cortical 
and  perhaps  even  some  sub¬ 
cortical  changes  in  protein 
concentration  and  myoinositol. 
In  Figure  2B,  the  green  is  the 
MS  patient,  and  the  blue  is  the 

average  over  healthy  volunteers.  We  are  exceptionally  excited  by  this  initial  result  and  felt  it  important  to  share 
here  as  Task  5  will  indeed  prove  the  sensitivity  of  these  advanced  techniques  to  MS. 


-0.2  -0.1  0  0.1  0.2 
Asymmetry  at  2  ppm  (AU) 


-02  -0.1  0  0.1  0.2 
Asymmetry  at  3  ppm  (AU) 


-02  -0.1  0  0.1  0.2 
Asymmetry  at  3.5  ppm  (AU) 


Figure  2:  (A)  APT  maps  for  all  slices  derived  in  a  healthy  volunteer  and 
concomitant  test-retest  CEST  spectra  in  2  healthy  volunteers.  The  test- 
retest  is  over  the  whole  brain,  and  the  green/blue  and  magenta/black 
spectra  pairs  are  from  the  same  volunteers.  (B)  histogram  analysis  of 
segmented  GM  in  20  healthy  volunteers  and  1  MS  patient  with  clinically 
noted  cognitive  impairment.  Note  that  histograms  at  2ppm  and  3.5ppm 
show  a  downward  shift  of  the  MS  patient  relative  to  the  healthy  control 
indicating  initial  sensitivity  to  the  pathology  of  cortical  GM  damage. 


Glutamate  CEST  (GluCEST)  analysis  and  results 

GluCEST  analysis  proceeded  as  presented  in  (8,9).  We  performed  GluCEST  analysis  in  20  healthy  volunteers 
and  4  MS  patients  at  the  time  of  this  report,  though  only  one  had  been  analyzed  and  is  presented  here.  In  short, 
GluCEST-weighted  images  were  collected  for  a  single  slice  with  high  spectral  saturation  fidelity  (see  Task  2)  at 
a  slice  slightly  superior  to  the  corpus  callosum.  Sample  GluCEST  weighted  images  are  shown  as  a  function  of 
offset  frequency  in  Figure  3,  top  panel.  From  these  maps,  the  GluCEST  spectra  were  corrected  for  BO  and  B1 
in  the  fashion  presented  in  (9),  and  a  GluCEST  asymmetry  map  at  Aoj  =  3.0ppm  was  generated.  Figure  3, 
bottom  panel  shows  two  healthy  volunteers  and  one  patient  with  MS,  clinically  diagnosed  with  cognitive 
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impairment.  First,  it  can  be  seen  that  the 
GluCEST  maps  show  excellent  contrast 
between  WM  and  GM,  with  the  GM 
having  higher  GluCEST  signal  than  the 
WM  (expected).  What  is  exciting  to  note 
(and  in  general,  part  of  Task  5)  is  the 
visual  differences  between  the  patient 
and  the  two  controls  shown  here.  The 
MS  patient  shows  elevated  GluCEST 
signal  (black  arrow)  on  the  left  side,  but 
apparently  diminished  GluCEST  signal 
in  the  right  cortical  GM  (magenta  arrow) 
compared  to  healthy  subject  1  (magenta 
arrow).  This  seems  to  indicate,  at  least  at 
the  early  stages,  that  GluCEST  is 
detecting  cortical  GM  differences 
between  healthy  and  MS  patients. 
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Figure  3:  (top)  example  GluCEST-weighted  images  as  a 
function  of  offset  frequency  (bottom)  comparison  of  GluCEST 
maps  for  two  healthy  volunteers  and  one  MS  patient.  Note  the 
apparent  differences  between  the  MS  patient  and  healthy 
volunteers  (black  and  magenta  arrows) 


As  with  the  APT  CEST,  we  examined 
the  entire  20  healthy  control  cohort  in 
comparison  to  the  MS  case  and  Figure 
4A  shows  the  average  GluCEST  spectra 
derived  from  GM  and  WM  in  healthy 
patients  (blue  and  black,  respectively) 
compared  to  the  segmented  GM  in  the 
MS  case  (red).  It  can  be  seen  that  the 
spectral  quality  is  high  and  there  is  visual 
difference  between  the  spectra  for  healthy 
and  MS  GM.  Further,  we  analyzed  the 
histogram  of  GluCEST  signals  for  all  GM 

voxels  in  all  healthy  volunteers  and  compared  that  to  the  1  MS  patient  clinically  diagnosed  with  cognitive 
challenges  (Figure  4B).  As  with  the  APT  CEST,  it  can  be  seen  that  the  MS  patient  shows  a  downward  trend 
compared  to  the  healthy  volunteer  indicating  the  possibility  of  being  sensitive  to  cortical  GM  pathology,  which 

we  will  study  in  detail  in  Task 
5.  When  we  started  this 
project,  we  decided  against  an 
“all-in-one”  CEST  acquisition 
scheme  and  rather  have 
deployed  two  CEST 
acquisitions.  One  sensitive  to 
APT  (and  apparently 
myoinositol)  and  one  sensitive 
to  glutamate.  From  Figure 
4B,  and  in  comparison  to 
Figure  2B,  the  histograms 
from  the  APT-CEST  analysis 
at  Aw  =  3.0ppm  (glutamate) 
show  now  difference  between 
healthy  controls  and  the  1  MS 
patient,  however,  when 

utilizing  the  GluCEST  acquisition,  in  the  same  patient,  a  difference  can  be  appreciated.  The  rationale  for  this  is 
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Figure  4:  (A)  Average  GluCEST  spectra  for  WM  and  GM  in  healthy  (blue 
black)  and  MS  patient  (red)  GM.  (B)  Histogram  analysis  shows  that  the 
MS  has  a  downward  shift  of  the  GluCEST  signal  compared  to  the  healthy 
volunteer  (blue) 


that  the  sensitivity  to  exchanging  species  is  determined  by  the  power  of  the  RF  CEST  saturation.  For  the  APT, 
the  exchange  rate  is  on  the  order  of  20-100Hz,  whereas  for  glutamate  amines,  the  exchange  rate  is  faster  (50- 
200Hz).  Thus,  to  be  maximally  sensitive  to  both,  two  separate  pulse  powers  are  necessary.  We  discovered  this 
as  part  of  Task  2  in  the  phantoms  and  are  proud  to  note  that  it  was  the  right  choice  going  forward. 

Quantitative  Magnetization  Transfer  (qMT)  analysis  and  results 

High-resolution  selective  inversion  recovery  (SIR)  qMT  was  performed  in  the  same  cohort  as  for  CEST  and 
analyzed  according  to  (10)  and  given  in  Appendix  3.  In  short,  an  inversion  recovery  MRI  sequence  was 
performed  using  a  modified  inversion  pulse  that  is  relatively  insensitive  to  B1  and  BO  inhomogeneities.  The 
inversion  times  were  selected  to  sample  the  bi-exponential  recovery  known  to  exist  when  magnetization  transfer 
is  present.  For  every  voxel,  the  SIR  signal  equation  was  fit  to  the  recovery  curve  and  the  exchange  rate  (kmf), 
pool  size  ratio  (PSR)  and  longitudinal  relaxation  time  (Rif)  was  fit.  Appendix  3  provides  the  manuscript  that 
contains  details  of  the  pulse  sequence,  and  fitting  method. 


Here  we  report  the  initial  analyses  and  results  from  the  newly  deployed  method.  Figure  5,  left  panel  shows  the 
PSR,  Rif,  and  kmf  maps  from  a  representative  healthy  volunteer  and  one  patient  with  MS  and  concomitant 
cognitive  impairment.  From  the  PSR  maps,  it  can  be  seen  the  high  level  of  discrimination  between  WM 
(yellow/orange)  and  GM  (blue)  which  agrees  well  with  the  Rif  maps.  Two  things  should  be  noted  for  the  MS 
patient.  First,  the  WM  shows  a  globally  decreased  PSR  which  is  indicative  of  demyelination  across  the  entire 
slice,  while  the  Rif  and  kmf  maps  do  not  show  a  similar  pattern  (discussed  next).  Secondly,  when  looking  at 
the  GM,  it  is  not  apparent  that  the  patient  and  the  control  have  different  PSR  values,  yet  MS  patients  are  known 


Figure  5:  (left)  qMT-derived  maps  for  a  healthy  control  (top)  and  MS  patient  (bottom),  (right) 
segmentation  results  (top)  and  histogram  analyses  for  the  average  healthy  controls  (red)  and  the  MS  patient 
(blue).  Note  the  downward  shift  of  the  PSR  for  both  WM  and  GM  in  this  patient,  while  R1  and  kmf  are 
indistinguishable. 


to  have  myelin  loss  in  the  GM.  However,  Figure  5,  right  panel,  shows  a  histogram  analysis  of  over  the  healthy 
volunteers  (N  =  20)  and  1  MS  patient.  The  top  row  shows  the  segmented  WM  and  GM  for  an  example  healthy 
volunteer,  and  the  bottom  panels  show  the  histograms  for  WM  and  GM  for  all  of  the  qMT-derived  indices.  It 
can  be  seen  again  that  for  WM,  the  MS  patient  has  a  substantially  downshifted  PSR,  normal  Rif  and  kmf. 

Importantly,  however,  in  the  GM,  the  patient  also  shows  a  small  downshift  of  the  PSR,  while  Rif  and  kmf  are 
not  markedly  different  (it  should  be  noted  that  kmf  showed  some  instabilities  in  this  patient).  This  is  important 


to  note  in  that  one  argument  about  MT  imaging  is  that  it  is  hypothesized  that  R1  drives  the  change  in  the  MT 
effect  moreso  than  does  the  macromolecular  content.  This  figure  shows  that  rather,  in  both  WM  and  GM,  the 
PSR  is  abnormal,  but  neither  R1  nor  kmf  are  indicating  the  sensitivity  to  WM  and  GM  macromolecular 
pathology.  Task  5  will  explore  this  further  when  we  examine  a  larger  MS  cohort. 
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fMRI  analysis  and  results 

In  addition  to  the  quantitative  measures  that  have  already  been  shown  and  at  the  advice  of  our  mentor,  Dr.  Paul 
Newhouse,  a  neurocognition  expert,  we  added  3  fMRI  acquisitions  to  the  MRI  paradigm.  Those  three  methods 
were  an  N-Back  task,  a  Trail-making  task,  and  a  resting-state  fMRI  acquisition.  The  N-back  and  Trail-making 
tasks  are  important  as  they  are  also  performed  outside  the  MRI  scanner,  so  we  will  be  able  to  provide  direct 
correlations  between  what  is  performed  in  the  MRI  and  outside  of  the  MRI.  This  further  allows  us  to  directly, 

|  .  ,  and  non- invasively  probe  cognitive 

N-back  3  Trail  MakinQ  performance  in  a  manner  that  is  not  only 

unique,  but  it  has  not  been  performed  in  the 
MS  population  at  IT.  We  are  encouraged  by 
the  initial  results  and  wish  to  present  those 
here.  We  performed  the  fMRI  in  20  healthy 
volunteers  and  1  MS  patient  at  the  time  of 
this  report.  Figure  6  shows  a  direct 
comparison  of  the  N-back  3  (left  panels)  and 
the  Trail-making  (right  panels)  for  a  single 
healthy  volunteer  (top  panels)  and  a  patient 
with  MS  (bottom  panels).  As  it  can  be  seen, 
at  the  same  significance  threshold,  there  are 

_  _ _ activation  differences  between  the  healthy 

Healthy  volunteer  and  the  MS  patient.  These  are 

- -  especially  noted  for  areas  in  the  superior 

.  _  cortex  where  working  memory  is  targeted.  A 
greater  confidence  will  be  gained  with  a 
group  analysis,  but  this  will  be  reserved  for 
Task  5. 

Figure  6  (left)  N-back  3  and  (right) 
Trail-making  fMRI  activation  patterns 
in  a  healthy  (top  panels)  volunteer  and 
one  patient  with  MS  (bottom  panels). 
Note  for  the  N-back,  there  is  less 
activation  in  the  MS  patient,  where  as 
for  the  Trail-making  task,  there  are 
greater  areas  of  activation. 
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Discussion  of  Task  4 

We  are  pleased  with  the  quality  of  data  that  has  been  generated  and  are  actively  enrolling  healthy  volunteer  sand 
patients  with  MS.  We  have  scanned  20  healthy  volunteers,  and  4  MS  patients  to  this  point,  and  have  had  8 
healthy  volunteers  return  for  a  2nd  visit.  We  are  slightly  behind  the  enrollment  expectation,  but  this  is  in  large 
part  due  to  the  ground  work  necessary  to  start  the  human  studies.  We  are,  however,  encouraged  that  because  of 
this  extensive  focus  on  sequence  optimization,  that  the  data  quality  remains  superb  in  all  subjects.  We  have 
scheduled  the  remainder  of  the  healthy  subjects,  developed  a  robust  pipeline  for  analysis,  and  have  shown  initial 
success  in  implementing  these  acquisitions  in  patients  with  MS.  We  expect  no  delays  in  finalizing  enrollment 
and  will  continue  with  reproducibility  analysis  from  which  we  can  gauge  the  level  of  expected  deviation  from 
normal  in  patients  with  MS. 
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Task  5  —  Implementation  in  Patients  with  MS  and  concomitant  cognitive  impairment 

The  objective  of  this  task  is  to  deploy,  and  analyze  the  MR1  acquisitions  in  patients  with  MS. 

We  have  implemented  the  MRI  paradigm  in  4  patients  with  MS  and  expect  no  delays  in  enrollment  for  the 
remainder  of  the  MS  patients.  The  details  of  preliminary  results  are  given  in  Task  4  with  the  goal  of  comparing 
both  image  quality,  and  quantitative  differences.  We  have  nothing  else  to  report  for  Task  5  at  the  time  of  this 
annual  report. 

Task  6  —  Cross-sectional  analysis  of  derived  indices  between  patients  with  MS  and  healthy  volunteers  and 
correlation  with  clinical  measures  of  cognitive  impairment  derived  from  the  Minimal  Assessment  of 
Cognitive  Function  in  MS  (MACFIMS) 

The  objective  of  this  task  is  to  compare  quantitative  MRI  indices  across  cohorts,  implement  neuropsychiatric 
evaluations  in  both  healthy  and  MS  patient  cohorts,  and  derive  correlations  with  quantitative  MR  indices. 

Neuropsychiatric  Assessment  Battery  (Outside  MRI  Scanner) 

In  collaboration  with  Dr.  Paul  Newhouse,  we  have  decided  to  additionally  obtain  neuropsychiatric  evaluations 
in  healthy  volunteers  in  addition  to  patients  with  MS.  This  will  provide  us  with  a  greater  understanding  of  the 
variance  across  cohorts.  Therefore,  we  have  obtained  neuropsychiatric  data  using  the  paradigm  below  in  20 
healthy  volunteers  (8  repeats)  and  4  MS  patients. 

Tasks  -  Outside  the  scanner  and  BEFORE  coming  to  Vanderbilt 

Questionnaire  and  survey  already  developed  in  REDCap  to  be  completed  at  home  and  in  a  calm  environment. 
These  surveys  will  collect  data  related  to  baseline  mood,  anxiety,  and  cognitive  profile. 

Tasks  -  Outside  the  scanner  at  Vanderbilt  (<  1  hour  total) 

•  Short  measure  of  day-of  mood/anxiety 

•  N-back  test  (2-back  or  3-back):  measures  working  memory 

•  PAS  AT:  measures  working  memory 

•  Trail  making  test  (both  A  and  B):  measures  planning/executive  function 

•  "Black  Box”  (choice  reaction  time,  critical  flicker  fusion;  pre-scan  and  post-scan):  measures  processing 
speed/reaction  time 

•  *Buschke  selective  reminding  test  (8  trials):  measures  include  encoding  and  long-term  memory 

•  Digit  Symbol  Substitution  Test/DDST:  measures  visual  memory 

•  Posner  cueing  task:  measures  attentional  shift 

KEY  RESEARCH  ACCOMPLISHMENTS: 

1 .  Developed  a  Chemical  Exchange  Saturation  Transfer  (CEST)  simulation  pipeline  to  model  the  effects  of 
glutamate,  amide  proton  transfer,  and  myoinositol  in  the  gray  matter  (GM)  at  7T.  These  simulations 
reflect  contributions  from  the  metabolite  of  interest,  the  magnetization  transfer  (MT)  effect,  and  the  direct 
water  saturation. 

2.  Developed  a  quantitative  magnetization  transfer  (qMT)  simulation  pipeline  to  model  the  effects  of 
macromolecular  concentration  of  the  white  matter  (WM)  and  GM  in  healthy  tissue  and  tissue  impacted 
by  multiple  sclerosis  (MS).  This  pipeline  incorporates  modeling  of  the  semisolid  fraction  while  taking 
into  account  relaxation  times  (T1  and  T2)  changes  that  are  known  to  occur  in  each  cohort. 

3.  Phantoms  have  been  created  that  attempt  to  model  the  in  vivo  scenario.  That  is,  they  have  varied 
concentration  of  metabolite  and  concomitant  concentration  of  semi-solid  components  (MT)  at  constant 
pH.  These  phantoms  include:  glutamate,  myoinositol,  glycogen,  bovine  serum  albumin  (BSA),  and 
agarose.  These  phantoms  can  be  leveraged  for  greater  understanding  of  the  in  vivo  results 


9 


4.  Developed  and  optimized  a  set  of  novel  MRI  acquisition  strategies  to  study  the  CEST  effects  of 
glutamate,  myoinositol,  amide  protons,  taking  into  account  corrections  for  B1  and  BO  inhomogeneities. 
These  have  been  deployed  in  vivo. 

5.  Developed  and  optimized  a  high-resolution  (1mm2  in-plane  resolution)  qMT  acquisition  that  is  of 
sufficient  resolution  to  assess  cortical  GM  and  juxtacortical  WM  in  healthy  and  MS  cohorts. 

6.  Developed  a  multi-modal  motion-correction,  coregistration  and  WM/GM/CSF  segmentation  strategy 
that  not  only  maps  the  confidence  of  the  measurements  made  in  small  structures,  but  also  puts  all 
acquisitions  in  the  same  space  for  descriptive  statistics  on  each  cohort. 

7.  Implemented  a  set  of  fMRI  experiments  to  assess  working  memory  and  resting-state  fluctuations  in 
patients  as  compared  to  healthy  controls.  This  will  work  in  conjunction  with  item  #8:  neuropsychiatric 
evaluation 

8.  Created  a  detailed  neuropsychiatric  evaluation  paradigm  for  assessing  cognitive  impairment  that  can  be 
performed  outside  of  the  scanner  environment  and  can  be  leveraged  for  correlation  testing  in  the  final 
year  of  this  project. 

9.  In  summary,  we  have  developed,  optimized,  and  deployed  a  multi-parametric,  multi-modal  MRI  toolkit 
to  assess  neurochemical,  macromolecular,  functional,  and  structural  changes  in  vivo  at  7T.  Additionally, 
we  have  developed  a  detailed  neuropsychiatric  evaluation  paradigm  to  utilize  for  comparisons  across 
cohorts,  which  can  be  further  extended  to  any  MRI  study  of  cognitive  impairment. 

CONCLUSION: 

We  have  concluded  the  first  year  of  this  project  and  have  many  significant  contributions  to  report.  First,  we  have 
for  the  first  time,  developed  a  battery  of  quantitative  MRI  methods  that  are  of  sufficient  resolution  and  sensitivity 
to  characterize  cortical  gray  matter  in  healthy  volunteers  and  patients  with  multiple  sclerosis.  To  this  end,  we 
have  a  <  1  hour  exam  card  that  can  be  deployed  on  any  7T  scanner  that  can  investigate  neurochemical 
composition,  macromolecular/myelin  deficiency,  and  functional  impairment  in  a  patient  cohort.  Additionally, 
we  have  generated  data  that  suggests  that  there  are  differences  between  healthy  volunteers  and  patients  with 
MS,  while  expected,  has  yet  to  be  shown  because  lower  MRI  field  strengths  have  insufficient  sensitivities  to 
these  macromolecules  and  neurochemicals,  and  insufficient  resolution  to  study  only  gray  matter.  We  have 
additionally  partnered  with  a  neurocognitive  expert,  and,  with  his  help,  developed  a  novel  neuropsychiatric 
battery  to  assess  cognition  in  MS.  We  understand  that  these  techniques,  while  not  novel,  have  not  been 
implemented  in  patients  with  MS,  and  may  provide  evidence  for  greater  scope  in  any  patient  with  neurocognitive 
decline.  We  have  studied  20  healthy  volunteers  and  a  handful  of  MS  patients  to  this  point  and  will  expand  the 
patient  and  control  enrollment  in  year  2.  We  have  finally,  developed  a  pipeline  for  analysis  that  requires 
minimal  human  interaction  and  will  deploy  this  for  real-time  analysis  in  year  2.  It  is  important  to  note,  however, 
that  while  these  techniques  are  currently  being  explored  for  use  at  7T,  in  year  2,  we  will  develop  a  lower  field, 
and  thus,  significantly  more  clinically  relevant,  set  of  exams  that  will  provide  a  similar  MRI  toolkit.  Lastly,  the 
exam  as  developed  here  is  not  specific  for  MS  and  can  be  implemented  in  a  wide  range  of  patients  and 
volunteers  to  explore  the  neurochemistry,  functional  processing,  and  macromolecular  composition  of  cortical 
gray  matter. 

PUBLICATIONS,  ABSTRACTS,  AND  PRESENTATIONS: 

1 .  Lay  Press:  Nothing  to  report 

2.  Peer-Reviewed  Scientific  Journals:  Nothing  to  report 

3 .  Invited  Articles :  N othing  to  report 

4.  Abstracts:  7  abstracts  to  the  International  Society  for  Magnetic  Resonance  in  Medicine  (ISMRM) 
annual  conference  have  been  prepared  and  will  be  submitted  in  Q1Y2  (November  12,  2014  deadline) 


INVENTIONS,  PATENTS  AND  LICENSES: 
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Nothing  to  Report 


REPORTABLE  OUTCOMES: 

1 .  High-resolution,  optimized  7T  MRI  acquisition  strategies  (so-called  Exam  Card  in  Philips  language)  to 
quantitatively  evaluate  the  macromolecular,  metabolic,  functional,  and  structural  characteristics  that  can 
(and  currently  is)  be  implemented  in  healthy  controls  and  any  patient  cohort.  A  complete  listing  of  the 
MRI  acquisition  paradigm  (Exam  Card)  is  given  in  Appendix  1 . 

2.  Assessment  of  the  reproducibility  and  stability  of  each  measurement  over  time  is  currently  ongoing. 

3.  A  complete  set  of  neuropsychiatric  assessments,  some  of  which  are  completely  novel  in  patients  with 
MS 

4.  Analysis  pipeline  for  generation  of  quantitative  MRI-derived  indices.  The  pipeline  includes  motion 
correction,  multi-modal  image  co-registration,  and  WM/GM/CSF  segmentation  along  with  generation  of 
quantitative  indices  for  further  statistical  comparisons. 

5.  A  CEST  simulation  GUI  for  further  studies  of  the  CEST  effect  in  vivo. 

6.  Collection  of  experiments,  simulations,  and  phantom  studies  that  have  provided  evidence  for  the 
minimal  resolution  attainable  while  maximizing  sensitivity  to  change  in  patient  populations. 

7.  Complete  set  of  MRI  and  neuropsychiatric  data  in  20  healthy  volunteers  (with  8  additional  repeat  visits) 
and  4  patients  with  MS  and  clinically  confirmed  cognitive  impairment.  These  data  are  summarized  in 
the  Overall  Project  Summary. 

OTHER  ACHIEVEMENTS: 

1 .  Because  of  the  nature  of  the  experiments  performed  as  a  part  of  this  grant,  that  is  to  implement  the 
highest  resolution  quantitative  MRI  at  7T  in  WM  and  GM  of  healthy  participants  and  MS  patients,  we 
have  been  able  to  extend  these  tools  to  the  spinal  cord,  which  has  significant  impact  and  scope  for 
patients  with  other  forms  of  neurological  injury.  One  manuscript  ( Dula  AN,  Pawate  S,  Dethrage  LM, 
Conrad  BN,  Barry  RL,  Smith  SA.  CEST  of  the  Cervical  Spinal  Cord  at  7  Tesla.  Submitted  to  NMR  in 
Biomed  on  30-Sept-2014 )  has  already  been  submitted  on  the  results  from  this  extension  to  other  parts  of 
the  body. 

2.  The  phantoms  and  simulations  that  were  created  are  not  MRI  field  strength  dependent.  Therefore,  we 
have  changed  the  simulations  to  study  the  impact  of  transitioning  to  a  lower-field  strength  (i.e.  3T)  for 
greater  clinical  implementation.  We  have  begun  to  study  the  sensitivity  for  3T  utilization.  This  is 
critical,  as  it  was  noted  in  our  initial  application,  that  7T  MRI  scanners  are  not  directly  clinically 
impactful.  Thus,  we  have  been  able  to  create  an  MRI  acquisition  strategy  at  lower  field  strength  which 
will  be  deployed  in  year  2.  We  are  excited  about  the  possibility  of  reaching  a  greater  clinical 
community  with  the  studies  in  year  2  at  3T. 

3.  From  the  preliminary  results  generated  in  the  first  year  of  this  award,  Dr.  Pawate  (co-investigator)  is 
preparing  a  grant  submission  to  the  National  Multiple  Sclerosis  Society  in  February  2015  to  study 
cerebral  changes  in  primary  progressive  MS  (PPMS)  patients. 


REFERENCES:  List  all  references  pertinent  to  the  report  using  a  standard  journal  format  (i.e.  format  used  in 
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E  CEST_Reddy_GluCEST  1 1 :37.0 
E  Bl_Reddy_multiAngle  01:42.0 
E  B0_Reddy_mulitecho  00:03.9 
EqMT  High  Res  02:38.7 
E  FMRI  RESTI NGSTATE  08:34.0 
E  FMRI_TRAI  LMAKI NG  04:14.0 
E  FMRI_nback  08:30.0 
ETl_3D_TFE_quantGeo  00:55.7 
ET2star_multiEcho  01:02.7 

E  {B1860723-4F8F-476e-8075-D42C65706693}  (0) 
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Q  Hospital  (2)  |  O20141021  CESTfMRI  (13)  52:37.4  |  E SCOUT  SHC32  00:28.7 


1 NFO  PAGE 

GEOMETRY 

CONTRAST 

Total  scan  duration 

00:28.7 

Multi-transmit 

no 

Scan  type 

Imaging 

Rel.  signal  level  (%) 

100 

Nucleus 

HI 

Scan  mode 

M2D 

Act.  TR/TE  (ms) 

7.8/  4.9 

Coil  selection  1 

RX-lntf-1 

technique 

FFE 

ACQ  matrix  M  x  P 

256  x  128 

Xmit  Coil  selection 

MTX-Volume- 

+  ZOOM 

no 

ACQ  voxel  MPS  (mm) 

0.98/  1.95  / 

T/R 

Contrast  enhancement 

Tl 

10.0 

User  def  elem  sel 

no 

Acquisition  mode 

cartesian 

REC  voxel  MPS  (mm) 

0.98/  0.98/ 

element  selection 

All 

Fast  Imaging  mode 

TFE 

10.0 

connection 

conn-A 

shot  mode 

multishot 

Scan  percentage  (%) 

50 

Coil  selection  2 

RX-lntf-2 

TFE  factor 

64 

TFE  shots 

2 

element  selection 

All 

startup  echoes 

default 

TFE  dur.  shot  /  acq  (ms) 

1062.9/  501.9 

Dual  coil 

yes 

+TFE  followup  echoes 

0 

TFE  shot  interval  (ms) 

1063.369 

Multi  coil 

no 

shot  interval 

shortest 

Min.  Tl  delay 

287.8819 

CLEAR 

no 

profile  order 

linear 

Act.  WFS  (pix)  /  BW 
(Hz) 

3.513  /  288.4 

FOV  FH  (mm) 

250 

Echoes 

1 

AP  (mm) 

250 

partial  echo 

no 

Min.  WFS  (pix)  /  Max. 

BW  (Hz) 

stack  RL  (mm) 

50 

shifted  echo 

no 

Voxel  size  FH  (mm) 

0.9765625 

in-phase 

Min.  TR/TE  (ms) 

7.8/  2.8 

TE 

RF  avg  power  computed 
(W) 

2.744635 

AP  (mm) 

1.953125 

(ms) 

4.93426 

Slice  thickness  (mm) 

10 

Flip  angle  (deg) 

15 

SAR  /  head 

<  100  % 

Recon  voxel  size  (mm) 

0.9765625 

TR 

shortest 

Whole  body  /  level 

<0.1  W/kg/ 

Fold-over  suppression 

no 

Halfscan 

no 

normal 

Reconstruction  matrix 

256 

Water-fat  shift 

user  defined 

B1  rms 

1.32  uT 

SENSE 

no 

(pixels) 

3.5 

PNS  /  level  //  VUIIS  : 

35  %  /  normal 

k-t  BLAST 

no 

Shim 

default 

dortch  : 

Stacks 

3 

mDIXON 

no 

Sound  Pressure  Level 
(dB) 

26.85212 

current 

A 

Fat  suppression 

no 

type 

parallel 

no 

MOTIOr 

water  suppression 

M 

slices 

3 

TFE  prepulse 

invert 

Cardiac  synchronization 

no 

slice  gap 

user  defined 

slice  selection 

no 

Heart  rate  >  250  bpm 

no 

gap  (mm) 

10 

shared 

no 

Respiratory 

no 

slice  orientation 

sagittal 

delay 

user  defined 

compensation 

fold -over  direction 

AP 

(ms) 

800 

Navigator  respiratory 
comp 

no 

fat  shift  direction 

F 

PSIR 

no 

Flow  compensation 

no 

Slice  scan  order 

default 

+inv  pulse  type 

-Fdefault 

fMRI  echo  stabilisation 

no 

Stack  scan  order 

ascend 

MTC 

no 

Motion  smoothing 

no 

Move  table  per  stack 

no 

T2prep 

no 

NSA 

1 

Stack  alignment 

no 

Research  prepulse 

no 

DYN/AN< 

G 

Stack  display  order 

no 

Diffusion  mode 

no 

Angio  /  Contrast  enh. 

no 

PlanAlign 

no 

Elastography  mode 

no 

Quantitative  flow 

no 

REST  slabs 

0 

SAR  mode 

high 

Manual  start 

no 

Catheter  tracking 

no 

B1  mode 

default 

+Abuse  dynamic  loop 

no 

Interactive  positioning 

no 

SAR  Patient  data 

auto 

Dynamic  study 

no 

Allow  table  movement 

no 

PNS  mode 

low 

OFFC/AI 

MG 

Gradient  mode 

default 

Arterial  Spin  labeling 

no 

POST/  PRC 

DC 

Stacks 

3 

SofTone  mode 

no 

current 

A 

Preparation  phases 

auto 

A 

Stack  Offc.  AP 
(P=+mm) 

0 

Interactive  FO 

no 

SENSE  ref.  scan 

no 

RL  (L=+mm) 

o 

SmartPlan  survey 

no 

FH  (H  =+mm) 

o 

BO  field  map 

no 

Ang.  AP  (deg) 

o 

B1  field  map 

no 

RL  (deg) 

o 

MIP/MPR 

no 

FH  (deg) 

o 

Images 

M,  no,  no,  no 

Autoview  image 

M 

Calculated  images 

no,  no,  no,  no 

Reference  tissue 

Grey  matter 

Preset  window  contrast 

soft 

Reconstruction  mode 

real  time 

Save  raw  data 

no 

Hardcopy  protocol 

no 

Ringing  filtering 

rectangular 

Geometry  correction 

default 
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Q  Hospital  (2)  |  O20141021  CESTfMRI  (13)  52:37.4  |  EwiP  MTX  SENSE  32ch  01:28.9 


1 NFO  PAGE 

GEOMETRY 

CONTRAST 

Total  scan  duration 

01:28.9 

Multi-transmit 

no 

Scan  type 

Imaging 

Rel.  signal  level  (%) 

100 

Nucleus 

HI 

Scan  mode 

3D 

Act.  TR/TE  (ms) 

8.0/  0.75 

Coil  selection  1 

RX-lntf-1 

technique 

FFE 

ACQ  matrix  MxP 

96  x  75 

Xmit  Coil  selection 

MTX-Volume- 

T/R 

loop  order 

zy order 

ACQ  voxel  MPS  (mm) 

5.52  /  7.07  / 

6.00 

+  ZOOM 

no 

User  def  elem  sel 

no 

Contrast  enhancement 

T1 

REC  voxel  MPS  (mm) 

5.52  /  5.52  / 

3.00 

element  selection 

All 

Acquisition  mode 

cartesian 

connection 

conn-A 

Fast  Imaging  mode 

none 

Scan  percentage  (%) 

78.125 

Coil  selection  2 

RX-lntf-2 

3D  non-selective 

no 

Packages 

1 

element  selection 

All 

Echoes 

1 

Act.  WFS  (pix)  /  BW 
(Hz) 

0.489  /  2071.3 

Dual  coil 

yes 

partial  echo 

no 

Multi  coil 

yes 

shifted  echo 

no 

Min.  WFS  (pix)  /  Max. 

BW  (Hz) 

0.486  /  2083.3 

CLEAR 

no 

TE 

shortest 

FOV  FH  (mm) 

530 

Flip  angle  (deg) 

1 

RF  avg  power  computed 
(W) 

0.05085949 

RL  (mm) 

530 

TR 

shortest 

SAR  /  head 

<  2  % 

stack  AP  (mm) 

300 

Halfscan 

no 

Whole  body  /  level 

0.0  W/kg  / 
normal 

Voxel  size  FH  (mm) 

5.520833 

Water -fat  shift 

minimum 

RL  (mm) 

7.066667 

Shim 

default 

B1  rms 

0.18  uT 

AP  (mm) 

3 

mDIXON 

no 

PNS  /  level  //  VUIIS  : 
dortch  : 

18  %  /  normal 

Recon  voxel  size  (mm) 

5.520833 

Fat  suppression 

no 

Image  shutter 

yes 

Water  suppression 

no 

Sound  Pressure  Level 
(dB) 

27.92985 

Fold-over  suppression 

no 

MTC 

no 

Slice  oversampling 

default 

Research  prepulse 

no 

MOTI  or 

\l 

RF  select.  FOS 

no 

Diffusion  mode 

no 

Cardiac  synchronization 

no 

Reconstruction  matrix 

96 

Elastography  mode 

no 

Heart  rate  >  250  bpm 

no 

SENSE 

no 

SAR  mode 

high 

Respiratory 

compensation 

no 

k-t  BLAST 

no 

B1  mode 

default 

Overcontiguous  slices 

yes 

SAR  Patient  data 

auto 

Navigator  respiratory 
comp 

no 

Stacks 

2 

PNS  mode 

low 

current 

A  slices 

Gradient  mode 

default 

Flow  compensation 

no 

100 

SofTone  mode 

no 

fMRI  echo  stabilisation 

no 

slice  orientation 

coronal 

NSA 

3 

SMART 

yes 

fold -over  direction 

RL 

fat  shift  direction 

F 

DYN/ ANG 

Chunks 

1 

Angio  /  Contrast  enh. 

no 

Stacks  as  packages 

no 

Quantitative  flow 

no 

Manual  start 

no 

Move  table  per  stack 

no 

Stack  alignment 

no 

+Abuse  dynamic  loop 

no 

Stack  display  order 

no 

Dynamic  study 

no 

PlanAlign 

no 

Arterial  Spin  labeling 

no 

REST  slabs 

0 

POST/  PROC 

Catheter  tracking 

no 

Preparation  phases 

full 

Interactive  positioning 

no 

1  nteractive  F0 

no 

Allow  table  movement 

no 

SENSE  ref.  scan 

yes 

OFFC/ANG 

SmartPlan  survey 

no 

Stacks 

2 

BO  field  map 

no 

current 

A 

B1  field  map 

no 

Stack  Offc.  AP 
(P=+mm) 

11.5916 

MIP/MPR 

no 

Images 

no,  no,  no,  no 

RL  (L=+mm) 

-3.434925 

Autoview  image 

no 

FH  (H  =+mm) 

4.208135 

Calculated  images 

no,  no,  no,  no 

Ang.  AP  (deg) 

0 

Reference  tissue 

White  matter 

RL  (deg) 

0 

Preset  window  contrast 

soft 

FH  (deg) 

0 

Reconstruction  mode 

immediate 

Save  raw  data 

no 

Hardcopy  protocol 

no 

Ringing  filtering 

default 

Elliptical  k-space  shutter 

default 
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Q  Hospital  (2)  |  O20141021  CESTfMRI  (13)  52:37.4  |  (Er  Tl_3D_TFE_isol.25mm_s2. 5s  SENSE_ Sagittal  02:12.2 


1 NFO  PAGE 

GEOMETRY 

CONTRAST 

Total  scan  duration 

02:12.2 

Multi-transmit 

no 

Scan  type 

Imaging 

Rel.  signal  level  (%) 

100 

Nucleus 

HI 

Scan  mode 

3D 

Act.  TR/TE  (ms) 

2.8/  1.32 

Coil  selection  1 

RX-lntf-1 

technique 

FFE 

ACQ  matrix  M  x  P 

204  x  204 

Xmit  Coil  selection 

MTX-Volume- 

+  ZOOM 

no 

ACQ  voxel  MPS  (mm) 

1.25  /  1.25  / 

T/R 

Contrast  enhancement 

Tl 

1.25 

User  def  elem  sel 

no 

Acquisition  mode 

cartesian 

REC  voxel  MPS  (mm) 

1.14/  1.14/ 

element  selection 

All 

Fast  Imaging  mode 

TFE 

1.25 

connection 

conn-A 

3D  non -selective 

no 

Scan  percentage  (%) 

100 

Coil  selection  2 

RX-lntf-2 

shot  mode 

multishot 

TFE  shots 

30 

element  selection 

All 

TFE  factor 

256 

TFE  dur.  shot  /  acq  (ms) 

1669.4/  716.5 

Dual  coil 

yes 

3D  free  factor 

no 

Min.  Tl  delay 

381.0045 

CLEAR 

yes 

startup  echoes 

default 

Act.  WFS  (pix)  /  BW 

0.579  /  1750.7 

body  tuned 

yes 

+TFE  followup  echoes 

0 

(Hz) 

FOV  FH  (mm) 

256 

shot  interval 

user  defined 

Min.  WFS  (pix)  /  Max. 

BW  (Hz) 

0.575  /  1763.2 

AP  (mm) 

256 

(ms) 

4500 

RF  avg  power  computed 
(W) 

0.8993402 

RL  (mm) 

172.5 

profile  order 

linear 

Voxel  size  FH  (mm) 

1.25 

turbo  direction 

radial 

SAR  /  head 

<  33  % 

AP  (mm) 

1.254902 

CENTRA  (spiral) 

no 

Whole  body  /  level 

0.0  W/kg  / 

RL  (mm) 

1.25 

Echoes 

1 

normal 

Recon  voxel  size  (mm) 

1.142857 

partial  echo 

no 

B1  rms 

0.76  uT 

Fold-over  suppression 

no 

shifted  echo 

no 

PNS  /  level  //  VUIIS  : 

57  %  /  normal 

Slice  oversampling 

default 

TE 

shortest 

dortch  : 

RF  select.  FOS 

no 

Flip  angle  (deg) 

7 

Sound  Pressure  Level 
(dB) 

37.43258 

Reconstruction  matrix 

224 

TR 

shortest 

SENSE 

yes 

Ha  If  scan 

no 

MOTI  or 

\l 

P  reduction  (AP) 

2 

Water -fat  shift 

minimum 

Cardiac  synchronization 

no 

P  os  factor 

1 

Shim 

auto 

Heart  rate  >  250  bpm 

no 

S  reduction  (RL) 

2 

mDIXON 

no 

Respiratory 

compensation 

no 

k-t  BLAST 

no 

Fat  suppression 

no 

Overcontiguous  slices 

no 

Water  suppression 

no 

Navigator  respiratory 
comp 

Stacks 

1 

TFE  prepulse 

invert 

Flow  compensation 

no 

slices 

138 

slice  selection 

no 

fMRI  echo  stabilisation 

no 

slice  orientation 

sagittal 

delay 

user  defined 

Motion  smoothing 

no 

fold -over  direction 

AP 

(ms) 

1300 

NSA 

1 

fat  shift  direction 

F 

PSIR 

no 

DYN/  ANi 

G 

Chunks 

1 

+inv  pulse  type 

+B1  opt  (low 
BW) 

PlanAlign 

no 

Angio  /  Contrast  enh. 

no 

REST  slabs 

0 

MTC 

no 

Quantitative  flow 

no 

no 

T2prep 

no 

CENTRA 

no 

Latneter  tracking 

Interactive  positioning 

no 

Research  prepulse 

no 

Manual  start 

no 

Allow  table  movement 

no 

Diffusion  mode 

no 

+Abuse  dynamic  loop 

no 

Elastography  mode 

no 

Dynamic  study 

no 

urru  AN  Li 

Stacks 

1 

SAR  mode 

high 

Arterial  Spin  labeling 

no 

Stack  Offc.  AP 

11.5916 

POST/  PROC 

di  moae 

aeTdUit 

(P=+mm) 

SAR  Patient  data 

auto 

Preparation  phases 

auto 

RL  (L=+mm) 

-3.434925 

PNS  mode 

low 

no 

i  nteractive  fu 

FH  (H  =+mm) 

4.208135 

Gradient  mode 

full  control 

CCMCC  rr\f 

no 

btNbt  rGT.  scan 

Ang.  AP  (deg) 

o 

max  strength 

33 

SmartPlan  survey 

no 

RL  (deg) 

0 

(mT/m) 

BO  field  map 

no 

FH  (deg) 

o 

max  slew  rate 

166 

B1  field  map 

no 

(T/m/s) 

MIP/MPR 

no 

Images 

M,  no,  no,  no 

Autoview  image 

M 

Calculated  images 

no,  no,  no,  no 

Reference  tissue 

Grey  matter 

Preset  window  contrast 

soft 

Reconstruction  mode 

immediate 

Save  raw  data 

no 

Hardcopy  protocol 

no 

Ringing  filtering 

rectangular 

Geometry  correction 

default 

Elliptical  k-space  shutter 

default 
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1 NFO  PAGE 

GEOMETRY 

CONTRAST 

Total  scan  duration 

09:09.6 

Multi-transmit 

no 

Scan  type 

Imaging 

Rel.  signal  level  (%) 

100 

Nucleus 

HI 

Scan  mode 

3D 

Act.  TR/TE  (ms) 

65  /  7.2 

Coil  selection  1 

RX-lntf-1 

technique 

FFE 

Dyn.  scan  time 

00:08.579 

Xmit  Coil  selection 

MTX-Volume- 

T/R 

+  ZOOM 

no 

Time  to  kO 

00:04.7 

Contrast 

enhancement 

T1 

ACQ  matrix  MxP 

160  x  148 

User  def  elem  sel 

no 

ACQ  voxel  MPS  (mm) 

1.50/  1.62  / 

10.0 

element  selection 

All 

Acquisition 

mode 

cartesian 

connection 

conn-A 

REC  voxel  MPS  (mm) 

0.94/  0.94/ 

5.00 

Coil  selection  2 

RX-lntf-2 

Fast  Imaging 
mode 

EPI 

element  selection 

All 

3D  non- 
selective 

no 

Scan  percentage  (%) 

92.77109 

Dual  coil 

yes 

Act.  WFS  (pix)  /  BW 
(Hz) 

8.715/  116.3 

CLEAR 

yes 

shot  mode 

multishot 

body  tuned 

no 

EPI  factor 

7 

BW  in  EPI  freq.  dir.  (Hz) 

1230.8 

FOV  AP  (mm) 

240 

epi 

direction 

Y 

Min.  WFS  (pix)  /  Max. 

BW  (Hz) 

8.664  /  117.0 

RL  (mm) 

240 

Min.  TR/TE  (ms) 

54/  5.9 

FH  (mm) 

165 

Echoes 

1 

Voxel  size  AP  (mm) 

1.5 

partial 

echo 

no 

RF  avg  power  computed 
(W) 

1.351341 

RL  (mm) 

1.5 

SAR  /  head 

<  49  % 

FH  (mm) 

5 

shifted 

echo 

no 

Whole  body  /  level 

<  0.1  W/kg/ 
normal 

Recon  voxel  size  (mm) 

0.9375 

Fold-over  suppression 

no 

TE 

user  defined 

B1  rms 

0.93  uT 

Slice  oversampling 

default 

(ms) 

7.2 

PNS  /  level  //  VUIIS  : 
dortch  : 

50  %  /  normal 

RF  select.  FOS 

no 

Flip  angle  (deg) 

5 

Reconstruction  matrix 

256 

TR 

user  defined 

Sound  Pressure  Level 
(dB) 

25.24222 

SENSE 

yes 

(ms) 

65 

P  reduction  (RL) 

2 

Ha  If  scan 

no 

MOTIOr 

VI 

P  os  factor 

1 

Water -fat  shift 

minimum 

Cardiac  synchronization 

no 

S  reduction  (FH) 

2 

Shim 

volume 

Heart  rate  >  250  bpm 

no 

k-t  BLAST 

no 

ShimAlign 

no 

Respiratory 

compensation 

no 

Overcontiguous  slices 

yes 

mDIXON 

no 

Stacks 

1 

Fat  suppression 

ProSet 

Navigator  respiratory 
comp 

no 

slices 

33 

pulse  type 

1331 

slice  orientation 

transverse 

Water 

suppression 

no 

Flow  compensation 

no 

fold -over  direction 

RL 

fMRI  echo  stabilisation 

no 

MTC 

+pulsed  qMT/CEST 

NSA 

1 

fat  shift  direction 

L 

+Duration 

(ms) 

25 

nVM /  AMI 

r: 

Chunks 

1 

PlanAlign 

no 

Angio  /  Contrast  enh. 

no 

+B1  Mode 

constant 

REST  slabs 

0 

Quantitative  flow 

no 

+Amp 

3 

Catheter  tracking 

no 

Manual  start 

no 

+B1  Units 

max  amp.  (uT) 

Interactive  positioning 

no 

+Abuse  dynamic  loop 

no 

+Pulses/TR 

1 

Allow  table  movement 

no 

Dynamic  study 

individual 

+Offset 

Mode 

file 

dyn  scans 

64 

Ul-I-W  AIN 

iij 

Stacks 

1 

+Offset 

File 

G:/ patch/ rf_offsets.txt 

recon  multiplier 

1 

dyn  scan  times 

shortest 

Stack  Offc.  AP 
(P=+mm) 

-5.532147 

+Offset 

Units 

PPm 

FOV  time  mode 

default 

RL  (L=+mm) 

-3.434925 

dummy  scans 

0 

FH  (H  =+mm) 

35.45897 

+RF  Shape 

gauss  cest  mt 

immediate 

subtraction 

no 

Ang.  AP  (deg) 

0 

Research 

prepulse 

no 

fast  next  scan 

no 

RL  (deg) 

-5.457352 

FH  (deg) 

0 

Diffusion  mode 

no 

synch,  ext.  device 

no 

Shim  Size  AP  (mm) 

180.7495 

Elastography 

mode 

no 

MTC 

no 

dyn  stabilization 

no 

RL  (mm) 

134.1432 

FH  (mm) 

114.2866 

SAR  mode 

high 

prospect,  motion 

corr. 

no 

Offc.  AP  (P=+mm) 

-1.87074 

B1  mode 

default 

SAR  Patient 
data 

auto 

Keyhole 

no 

RL  (L=+mm) 

-3.660397 

Arterial  Spin  labeling 

no 

FH  (H  =+mm) 

33.86363 

PNS  mode 

low 

Ang.  AP  (deg) 

2.26406 

POST/  PROC 

«  i  ■  i  i 

default 

RL  (deg) 

-6.580403 

uraaient  mooe 

Preparation  phases 

sameprep 

CrvfTrvno  rv\ nrlo 

no 

Interactive  FO 

no 

FH  (deg) 

4.637685 

jot  i  one  mooe 

SENSE  ref.  scan 

no 

SmartPlan  survey 

no 

BO  field  map 

no 

B1  field  map 

no 

MIP/MPR 

no 

Images 

M,  no,  no,  no 

Autoview  image 

M 

Calculated  images 

no,  no,  no,  no 

Reference  tissue 

Grey  matter 

EPI  2D  phase  correction 

no 

Preset  window  contrast 

soft 

Reconstruction  mode 

real  time 

reuse  memory 

no 

Save  raw  data 

no 

Hardcopy  protocol 

no 

Ringing  filtering 

rectangular 

Geometry  correction 

default 

Elliptical  k-space  shutter 

default 
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1 NFO  PAGE 

GEOMETRY 

CONTRAST 

Total  scan  duration 

11:37.0 

Multi-transmit 

no 

Scan  type 

Imaging 

Rel.  signal  level  (%) 

100 

Nucleus 

HI 

Scan  mode 

2D 

Act.  TR/TE  (ms) 

5.6/  2.7 

Coil  selection  1 

RX-lntf-1 

technique 

FFE 

Dyn.  scan  time 

00:13.909 

Xmit  Coil  selection 

MTX-Volume- 

T/R 

+  ZOOM 

no 

Time  to  kO 

00:08.6 

Contrast  enhancement 

T1 

ACQ  matrix  MxP 

128  x 126 

User  def  elem  sel 

no 

Acquisition  mode 

cartesian 

ACQ  voxel  MPS  (mm) 

1.88/  1.90/ 

10.0 

element  selection 

All 

Fast  Imaging  mode 

TFE 

connection 

conn-A 

shot  mode 

multishot 

REC  voxel  MPS  (mm) 

0.94/  0.94/ 

10.0 

Coil  selection  2 

RX-lntf-2 

TFE  factor 

3 

element  selection 

All 

startup  echoes 

default 

Scan  percentage  (%) 

98.4375 

Dual  coil 

yes 

+TFE  followup  echoes 

0 

TFE  shots 

42 

CLEAR 

yes 

shot  interval 

shortest 

TFE  dur.  shot  /  acq  (ms) 

134.0/  16.8 

body  tuned 

no 

profile  order 

low high 

TFE  shot  interval  (ms) 

331.1945 

FOV  AP  (mm) 

240 

Echoes 

1 

Act.  WFS  (pix)  /  BW 
(Hz) 

0.700  /  1446.8 

RL  (mm) 

240 

partial  echo 

no 

Voxel  size  AP  (mm) 

1.87 

shifted  echo 

no 

Min.  WFS  (pix)  /  Max. 

BW  (Hz) 

0.649  /  1562.5 

RL  (mm) 

1.87 

TE 

user  defined 

Min.  TR/TE  (ms) 

4.0/  1.36 

Slice  thickness  (mm) 

10 

(ms) 

2.7 

Recon  voxel  size  (mm) 

0.9375 

Flip  angle  (deg) 

10 

RF  avg  power  computed 
(W) 

2.058919 

Fold-over  suppression 

no 

TR 

user  defined 

SAR  /  head 

<  75  % 

Reconstruction  matrix 

256 

(ms) 

5.6 

Whole  body  /  level 

<  0.1  W/kg/ 
normal 

SENSE 

no 

Ha  If  scan 

no 

k-t  BLAST 

no 

Water -fat  shift 

minimum 

B1  rms 

1.15  uT 

Slice  orientation 

transverse 

Shim 

volume 

PNS  /  level  //  VUIIS  : 
dortch  : 

31  %  /  normal 

Fold-over  direction 

RL 

ShimAlign 

no 

Fat  shift  direction 

P 

mDIXON 

no 

Sound  Pressure  Level 
(dB) 

17.13024 

PlanAlign 

no 

Fat  suppression 

no 

REST  slabs 

0 

Water  suppression 

no 

MOTI  ON 

Catheter  tracking 

no 

TFE  prepulse 

no 

Cardiac  synchronization 

no 

Interactive  positioning 

no 

MTC 

+pulsed 

qMT/CEST 

Heart  rate  >  250  bpm 

no 

Allow  table  movement 

no 

Respiratory 

compensation 

no 

OFFC/AN 

IG 

+Duration  (ms) 

10 

Slice  Offc.  AP 
(P=+mm) 

-4.247866 

+B1  Mode 

constant 

Navigator  respiratory 
comp 

no 

+Amp 

4.25 

RL  (L=+mm) 

-3.434925 

+B1  Units 

max  amp.  (uT) 

Flow  compensation 

no 

FH  (H  =+mm) 

35.88707 

+Pulses/TR 

10 

fMRI  echo  stabilisation 

no 

Ang.  AP  (deg) 

0 

+Duty  Cycle 

0.9 

Motion  smoothing 

no 

RL  (deg) 

-5.457352 

+ Offset  Mode 

baseline+range 

NSA 

1 

FH  (deg) 

0 

+  Min  Offset 

-5 

DYN/ ANG 

Shim  Size  AP  (mm) 

180.7495 

+Max  Offset 

5 

Angio  /  Contrast  enh. 

no 

RL  (mm) 

134.1432 

+Offset  Units 

ppm 

Quantitative  flow 

no 

FH  (mm) 

106.6376 

+RF  Shape 

gauss cest mt 

Manual  start 

no 

Offc.  AP  (P=+mm) 

-1.87074 

+ Interpulse 
Spoiling 

no 

+Abuse  dynamic  loop 

no 

RL  (L=+mm) 

-3.660397 

Dynamic  study 

individual 

FH  (H  =+mm) 

33.86363 

T2prep 

no 

dyn  scans 

50 

Ang.  AP  (deg) 

2.26406 

Research  prepulse 

no 

recon  multiplier 

1 

RL  (deg) 

-6.580403 

Diffusion  mode 

no 

dyn  scan  times 

shortest 

FH  (deg) 

4.637685 

Elastography  mode 

no 

FOV  time  mode 

default 

SAR  mode 

moderate 

dummy  scans 

0 

B1  mode 

default 

immediate 

subtraction 

no 

SAR  Patient  data 

auto 

fast  next  scan 

no 

PNS  mode 

low 

Gradient  mode 

default 

synch,  ext.  device 

no 

SofTone  mode 

no 

MTC 

no 

dyn  stabilization 

no 

prospect,  motion 

corr. 

no 

Keyhole 

no 

Arterial  Spin  labeling 

no 

POST/  PROC 

Preparation  phases 

sameprep 

Interactive  FO 

no 

SENSE  ref.  scan 

no 

SmartPlan  survey 

no 

BO  field  map 

no 

B1  field  map 

no 

MIP/MPR 

no 

Images 

M,  no,  no,  no 

Autoview  image 

M 

Calculated  images 

no,  no,  no,  no 

Reference  tissue 

Grey  matter 

Preset  window  contrast 

soft 

Reconstruction  mode 

real  time 

reuse  memory 

no 

Save  raw  data 

no 

Hardcopy  protocol 

no 

Ringing  filtering 

rectangular 

Geometry  correction 

default 
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1 NFO  PAGE 

GEOMETRY 

CONTRAST 

Total  scan  duration 

01:42.0 

Multi-transmit 

no 

Scan  type 

Imaging 

Rel.  signal  level  (%) 

100 

Nucleus 

HI 

Scan  mode 

MS 

Act.  TR  (ms) 

6000 

Coil  selection  1 

RX-lntf-1 

technique 

SE 

Act.  TE  (ms) 

12 

Xmit  Coil  selection 

MTX-Volume- 

T/R 

+  ZOOM 

no 

Dyn.  scan  time 

00:48.000 

Modified  SE 

no 

ACQ  matrix  MxP 

128  x 120 

User  def  elem  sel 

no 

+Optim  ref/crush 

no 

ACQ  voxel  MPS  (mm) 

1.88/  2.00/ 

10.0 

element  selection 

All 

Acquisition  mode 

cartesian 

connection 

conn-A 

Fast  Imaging  mode 

TSE 

REC  voxel  MPS  (mm) 

0.94/  0.94/ 

10.0 

Coil  selection  2 

RX-lntf-2 

shot  mode 

multishot 

element  selection 

All 

TSE  factor 

15 

Scan  percentage  (%) 

93.75 

Dual  coil 

yes 

startup  echoes 

0 

Packages 

1 

CLEAR 

no 

-Ffollowup  echoes 

0 

Min.  slice  gap  (mm) 

10 

FOV  AP  (mm) 

240 

profile  order 

lowhigh 

WFS(pix)/  BW  (Hz) 

1.986  /  510.1 

RL  (mm) 

240 

DRIVE 

no 

TSE  es  /  shot  (ms) 

12.0/  180 

FH  (mm) 

10 

ultrashort 

no 

TEeff  /  TEequiv  (ms) 

12  /  12 

Voxel  size  AP  (mm) 

1.87 

strong  FID 
crushing 

no 

Min.  TR  (ms) 

589 

RL  (mm) 

1.87 

RF  avg  power  computed 
(W) 

0.2697158 

Slice  thickness  (mm) 

10 

Echoes 

1 

Recon  voxel  size  (mm) 

0.9375 

partial  echo 

no 

SAR  /  head 

<  10  % 

Small  FOV  imaging 

no 

TE 

user  defined 

Whole  body  /  level 

0.0  W/kg  / 

Fold-over  suppression 

no 

(ms) 

12 

B1  rms 

0.41  uT 

Reconstruction  matrix 

256 

Flip  angle  (deg) 

60 

SENSE 

no 

Refocusing  control 

no 

PNS  /  level  //  VUIIS  : 
dortch  : 

30  %  /  normal 

k-t  BLAST 

no 

TR 

user  defined 

Sound  Pressure  Level 
(dB) 

16.02466 

Stacks 

1 

(ms) 

6000 

type 

parallel 

Ha  If  scan 

no 

+  Ref  Pulse  Shape 

slices 

1 

Water -fat  shift 

user  defined 

+  Ref  Pulse  Dur  [msec] 

slice  gap 

user  defined 

(pixels) 

2 

+  Ref  MAX  Dephase 

gap  (mm) 

0 

Shim 

volume 

+  Ref  MIN  Dephase 

slice  orientation 

transverse 

ShimAlign 

no 

+Crusher  b  value 

fold -over  direction 

RL 

mDIXON 

no 

MOTION 

fat  shift  direction 

P 

Fat  suppression 

no 

Cardiac  synchronization 

no 

Minimum  number  of 
packages 

1 

Water  suppression 

no 

Grad.  rev.  offres.  supp. 

no 

Heart  rate  >  250  bpm 

no 

Slice  scan  order 

default 

BB  pulse 

no 

Respiratory 

compensation 

no 

PlanAlign 

no 

MTC 

no 

Navigator  respiratory 
comp 

no 

REST  slabs 

0 

T2prep 

no 

Catheter  tracking 

no 

Research  prepulse 

no 

Flow  compensation 

no 

Interactive  positioning 

no 

Zoom  imaging 

no 

Temporal  slice  spacing 

default 

Allow  table  movement 

no 

Diffusion  mode 

no 

Motion  smoothing 

no 

OFFC/ANG 

Elastography  mode 

no 

NSA 

1 

Stacks 

1 

SAR  mode 

high 

DYN/ANG 

Stack  Offc.  AP 
(P=+mm) 

-4.247866 

B1  mode 

user  defined 

Manual  start 

no 

amplitude  (uT) 

9.5 

+Abuse  dynamic  loop 

flip  angle 
diminish 

RL  (L=+mm) 

-3.434925 

SAR  Patient  data 

auto 

FH  (H  =+mm) 

35.88707 

PNS  mode 

high 

Dynamic  study 

individual 

Ang.  AP  (deg) 

0 

Gradient  mode 

default 

dyn  scans 

2 

RL  (deg) 

-5.457352 

SofTone  mode 

no 

recon  multiplier 

1 

FH  (deg) 

0 

dyn  scan  times 

shortest 

Shim  Size  AP  (mm) 

180.7495 

FOV  time  mode 

default 

RL  (mm) 

134.1432 

dummy  scans 

0 

FH  (mm) 

106.6376 

immediate 

subtraction 

no 

Offc.  AP  (P=+mm) 

-1.87074 

RL  (L=+mm) 

-3.660397 

fast  next  scan 

no 

FH  (H  =+mm) 

33.86363 

synch,  ext.  device 

no 

Ang.  AP  (deg) 

2.26406 

dyn  stabilization 

no 

RL  (deg) 

-6.580403 

prospect,  motion 

corr. 

no 

FH  (deg) 

4.637685 

Keyhole 

no 

Arterial  Spin  labeling 

no 

POST/  PROC 

Preparation  phases 

sameprep 

Interactive  FO 

no 

SENSE  ref.  scan 

no 

SmartPlan  survey 

no 

BO  field  map 

no 

B1  field  map 

no 

MIP/MPR 

no 

Images 

M,  no,  no,  no 

Autoview  image 

M 

Calculated  images 

no,  no,  no,  no 

Reference  tissue 

Grey  matter 

Preset  window  contrast 

soft 

Reconstruction  mode 

real  time 

reuse  memory 

no 

Save  raw  data 

no 

Hardcopy  protocol 

no 

Ringing  filtering 

rectangular 

Geometry  correction 

default 
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1 NFO  PAGE 

GEOMETRY 

CONTRAST 

Total  scan  duration 

00:03.9 

Multi-transmit 

no 

Scan  type 

Imaging 

Rel.  signal  level  (%) 

100 

Nucleus 

HI 

Scan  mode 

2D 

Act.  TR/TEl/delta  TE 
(ms) 

53  /  3.4/  3.9 

Coil  selection  1 

RX-lntf-1 

technique 

FFE 

Xmit  Coil  selection 

MTX-Volume- 

T/R 

+  ZOOM 

no 

ACQ  matrix  M  x  P 

128  x 128 

Contrast  enhancement 

no 

ACQ  voxel  MPS  (mm) 

1.88/  1.88/ 

10.0 

User  def  elem  sel 

no 

Acquisition  mode 

cartesian 

element  selection 

All 

Fast  Imaging  mode 

none 

REC  voxel  MPS  (mm) 

0.94/  0.94/ 

10.0 

connection 

conn-A 

Echoes 

4 

Coil  selection  2 

RX-lntf-2 

partial  echo 

no 

Scan  percentage  (%) 

100 

element  selection 

All 

shifted  echo 

no 

Act.  WFS  (pix)  /  BW 
(Hz) 

2.824  /  358.8 

Dual  coil 

yes 

TE  first 

shortest 

Min.  WFS  (pix)  /  Max. 

BW  (Hz) 

0.649  /  1562.5 

CLEAR 

yes 

echospacing 

shortest 

body  tuned 

no 

flyback 

yes 

RF  avg  power  computed 
(W) 

2.058916 

FOV  AP  (mm) 

240 

Flip  angle  (deg) 

65 

RL  (mm) 

240 

TR 

shortest 

SAR  /  head 

<  75  % 

Voxel  size  AP  (mm) 

1.87 

Halfscan 

no 

Whole  body  /  level 

<0.1  W/kg/ 
normal 

RL  (mm) 

1.875 

Water -fat  shift 

maximum 

Slice  thickness  (mm) 

10 

Shim 

volume 

B1  rms 

1.15  uT 

Recon  voxel  size  (mm) 

0.9375 

ShimAlign 

no 

PNS  /  level  //  VUIIS  : 
dortch  : 

46  %  /  normal 

Fold-over  suppression 

no 

mDIXON 

no 

Reconstruction  matrix 

256 

Fat  suppression 

no 

Sound  Pressure  Level 
(dB) 

21.87152 

SENSE 

yes 

Water  suppression 

no 

P  reduction  (RL) 

2 

MTC 

no 

MOTION 

P  os  factor 

1 

Research  prepulse 

no 

Cardiac  synchronization 

no 

k-t  BLAST 

no 

Diffusion  mode 

no 

Heart  rate  >  250  bpm 

no 

Slice  orientation 

transverse 

Elastography  mode 

no 

Respiratory 

compensation 

no 

Fold-over  direction 

RL 

SAR  mode 

moderate 

Fat  shift  direction 

P 

B1  mode 

default 

Navigator  respiratory 
comp 

no 

PlanAlign 

no 

SAR  Patient  data 

auto 

Flow  compensation 

yes 

REST  slabs 

0 

PNS  mode 

low 

fMRI  echo  stabilisation 

no 

Catheter  tracking 

no 

Gradient  mode 

default 

NSA 

1 

Interactive  positioning 

no 

SofTone  mode 

no 

DYN/ANi 

G 

Allow  table  movement 

no 

Angio  /  Contrast  enh. 

no 

OFFC/ANG 

Slice  Offc.  AP 
(P=+mm) 

-4.247866 

Quantitative  flow 

no 

Manual  start 

no 

+Abuse  dynamic  loop 

no 

RL  (L=+mm) 

-3.434925 

FH  (H  =+mm) 

35.88707 

Dynamic  study 

no 

Ang.  AP  (deg) 

0 

Arterial  Spin  labeling 

no 

RL  (deg) 

-5.457352 

POST/  PROC 

FH  (deg) 

0 

Preparation  phases 

sameprep 

Shim  Size  AP  (mm) 

180.7495 

Interactive  FO 

no 

SENSE  ref.  scan 

no 

RL  (mm) 

134.1432 

FH  (mm) 

106.6376 

SmartPlan  survey 

no 

Offc.  AP  (P=+mm) 

-1.87074 

BO  field  map 

no 

RL  (L=+mm) 

-3.660397 

B1  field  map 

no 

FH  (H  =+mm) 

33.86363 

MIP/MPR 

no 

Ang.  AP  (deg) 

2.26406 

Images 

M,  R,  1 ,  no 

RL  (deg) 

-6.580403 

Autoview  image 

M 

FH  (deg) 

4.637685 

Calculated  images 

no,  no,  no,  no 

Reference  tissue 

Grey  matter 

Preset  window  contrast 

soft 

Reconstruction  mode 

real  time 

Save  raw  data 

yes 

Hardcopy  protocol 

no 

Ringing  filtering 

rectangular 

Geometry  correction 

default 
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Q  Hospital  (2)  |  O20141021  CESTfMRI  (13)  52:37.4  |  EqMTHigh  Res  02:38.7 


1 NFO  PAGE 

GEOMETRY 

CONTRAST 

Total  scan  duration 

02:38.7 

Multi-transmit 

no 

Scan  type 

Imaging 

Rel.  signal  level  (%) 

100 

Nucleus 

HI 

Scan  mode 

3D 

Act.  TR/TE  (ms) 

4.2/  2.2 

Coil  selection  1 

RX-lntf-1 

technique 

FFE 

Dyn.  scan  time 

00:11.251 

Xmit  Coil  selection 

MTX-Volume- 

T/R 

+  ZOOM 

no 

Time  to  kO 

00:07.0 

Contrast  enhancement 

T1 

ACQ  matrix  MxP 

212  x  210 

User  def  elem  sel 

no 

Acquisition  mode 

cartesian 

ACQ  voxel  MPS  (mm) 

1.00/  1.01/ 

2.00 

element  selection 

All 

Fast  Imaging  mode 

TFE 

connection 

conn-A 

3D  non -selective 

no 

REC  voxel  MPS  (mm) 

0.95  /  0.95  / 

2.00 

Coil  selection  2 

RX-lntf-2 

shot  mode 

multishot 

element  selection 

All 

TFE  factor 

54 

Scan  percentage  (%) 

99.08257 

Dual  coil 

yes 

3D  free  factor 

no 

TFE  shots 

12 

CLEAR 

yes 

startup  echoes 

user  defined 

TFE  dur.  shot  /  acq  (ms) 

233.8  /  225.8 

body  tuned 

no 

(number) 

0 

TFE  shot  interval  (ms) 

3353.226 

FOV  AP  (mm) 

212 

+TFE  followup  echoes 

0 

Act.  WFS  (pix)  /  BW 
(Hz) 

1.142  /  887.5 

RL  (mm) 

212 

shot  interval 

shortest 

FH  (mm) 

10 

profile  order 

low  high 

Min.  WFS  (pix)  /  Max. 

BW  (Hz) 

1.133  /  894.1 

Voxel  size  AP  (mm) 

1 

turbo  direction 

Y 

RF  avg  power  computed 
(W) 

1.347757 

RL  (mm) 

1 

Echoes 

1 

FH  (mm) 

2 

partial  echo 

no 

SAR  /  head 

<  49  % 

Recon  voxel  size  (mm) 

1 

shifted  echo 

no 

Whole  body  /  level 

<  0.1  W/kg/ 
normal 

Fold-over  suppression 

no 

TE 

shortest 

Slice  oversampling 

default 

Flip  angle  (deg) 

15 

B1  rms 

0.93  uT 

RF  select.  FOS 

no 

TR 

shortest 

PNS  /  level  //  VUIIS  : 
dortch  : 

50  %  /  normal 

Reconstruction  matrix 

224 

Halfscan 

no 

SENSE 

yes 

Water -fat  shift 

minimum 

Sound  Pressure  Level 
(dB) 

29.07384 

P  reduction  (RL) 

2 

Shim 

PB-volume 

P  os  factor 

1 

ShimAlign 

no 

MOTION 

S  reduction  (FH) 

1 

mDIXON 

no 

Cardiac  synchronization 

no 

k-t  BLAST 

no 

Fat  suppression 

no 

Heart  rate  >  250  bpm 

no 

Overcontiguous  slices 

no 

Water  suppression 

no 

Respiratory 

compensation 

no 

Stacks 

1 

TFE  prepulse 

no 

slices 

5 

MTC 

+SIR 

Navigator  respiratory 
comp 

no 

slice  orientation 

transverse 

+pulse  dur  (ms) 

5.5 

fold -over  direction 

RL 

+pulse  shape 

BRASORF 

Flow  compensation 

no 

fat  shift  direction 

P 

+ offset  (Hz) 

0 

fMRI  echo  stabilisation 

no 

Motion  smoothing 

no 

Chunks 

1 

-t-spoil  amp 
(mT/m) 

20 

PlanAlign 

no 

NSA 

1 

< 

I 

G 

REST  slabs 

0 

+spoil  dur  (ms) 

2 

Angio  /  Contrast  enh. 

no 

Catheter  tracking 

no 

+TFE  Saturation 

yes 

Interactive  positioning 

no 

+td  Mode 

constant 

Quantitative  flow 

no 

Allow  table  movement 

no 

+td  (ms) 

2500 

CENTRA 

no 

+ti  (ms) 

6,  10,  16,  26, 

42,  68,  110, 

178,  288,  468, 
760,  1233, 

2000,  8000, 

2000,  10000, 
1000,  1000, 

1000,  1000, 

1000,  1000, 

1000,  1000, 

1000,  1000, 

1000,  1000, 

1000,  1000, 

1000,  1000, 

1000,  1000, 

1000,  1000, 

1000,  1000, 

1000,  1000, 

1000,  1000, 

1000,  1000, 

1000,  1000, 

1000,  1000, 

1000,  1000, 

1000,  1000, 

1000,  1000, 

1000,  1000, 

1000,  1000, 

1000,  1000, 

1000,  1000, 

1000,  1000 

Manual  start 

no 

urrt/  aim 

l<J 

+Abuse  dynamic  loop 

no 

Stacks 

1 

Dynamic  study 

individual 

Stack  Offc.  AP 
(P=+mm) 

-4.247866 

dyn  scans 

14 

RL  (L=+mm) 

-3.434925 

recon  multiplier 

1 

FH  (H  =+mm) 

35.88707 

dyn  scan  times 

shortest 

Ang.  AP  (deg) 

0 

FOV  time  mode 

default 

RL  (deg) 

-5.457352 

dummy  scans 

0 

FH  (deg) 

0 

immediate 

subtraction 

no 

Shim  Size  AP  (mm) 

180.7495 

fast  next  scan 

no 

RL  (mm) 

134.1432 

synch,  ext.  device 

no 

FH  (mm) 

106.6376 

MTC 

no 

Offc.  AP  (P=+mm) 

-1.87074 

dyn  stabilization 

no 

RL  (L=+mm) 

-3.660397 

FH  (H  =+mm) 

33.86363 

prospect,  motion 

corr. 

no 

Ang.  AP  (deg) 

2.26406 

Keyhole 

no 

RL  (deg) 

-6.580403 

Arterial  Spin  labeling 

no 

FH  (deg) 

4.637685 

POST/  PR< 

DC 

Preparation  phases 

auto 

Interactive  FO 

no 

SENSE  ref.  scan 

no 

SmartPlan  survey 

no 

BO  field  map 

no 

B1  field  map 

no 

T2prep 

no 

MIP/MPR 

no 

Research  prepulse 

no 

Images 

M,  R,  1,  no 

Diffusion  mode 

no 

Autoview  image 

M 

Elastography  mode 

no 

Calculated  images 

no,  no,  no,  no 

SAR  mode 

high 

Reference  tissue 

Grey  matter 

B1  mode 

default 

Preset  window  contrast 

soft 

SAR  Patient  data 

auto 

Reconstruction  mode 

real  time 

PNS  mode 

high 

reuse  memory 

no 

Gradient  mode 

default 

Save  raw  data 

no 

SofTone  mode 

no 

Hardcopy  protocol 

no 

Ringing  filtering 

rectangular 

Geometry  correction 

default 

Elliptical  k-space  shutter 

default 
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Q  Hospital  (2)  |  O20141021  CESTfMRI  (13)  52:37.4  |  E  FMRI_RESTI  NGSTATE  08:34.0 


1 NFO  PAGE 

GEOMETRY 

CONTRAST 

Total  scan  duration 

08:34.0 

Multi-transmit 

no 

Scan  type 

Imaging 

Rel.  signal  level  (%) 

100 

Nucleus 

HI 

Scan  mode 

MS 

Act.  TR/TE  (ms) 

2000  /  25 

Coil  selection  1 

RX-lntf-1 

technique 

FFE 

Dyn.  scan  time 

00:02.000 

Xmit  Coil  selection 

MTX-Volume- 

T/R 

+  ZOOM 

no 

Time  to  kO 

00:15.0 

Contrast  enhancement 

no 

ACQ  matrix  M  x  P 

96x95 

User  def  elem  sel 

no 

Acquisition  mode 

cartesian 

ACQ  voxel  MPS  (mm) 

2.50/  2.50/ 

2.50 

element  selection 

All 

Fast  Imaging  mode 

EPI 

connection 

conn-A 

shot  mode 

single-shot 

REC  voxel  MPS  (mm) 

2.50  /  2.50  / 

2.50 

Coil  selection  2 

RX-lntf-2 

Echoes 

1 

element  selection 

All 

partial  echo 

no 

Scan  percentage  (%) 

100 

Dual  coil 

yes 

shifted  echo 

no 

Packages 

1 

CLEAR 

yes 

TE 

user  defined 

Min.  slice  gap  (mm) 

0 

body  tuned 

yes 

(ms) 

25 

EPI  factor 

37 

FOV  RL  (mm) 

240 

Flip  angle  (deg) 

63 

Act.  WFS  (pix)  /  BW 
(Hz) 

21.825  /  46.4 

AP  (mm) 

240 

TR 

user  defined 

FH  (mm) 

115 

(ms) 

2000 

BW  in  EPI  freq.  dir.  (Hz) 

2878.9 

Voxel  size  RL  (mm) 

2.5 

Ha  If  scan 

no 

Min.  WFS  (pix)  /  Max. 

BW  (Hz) 

21.770  /  46.5 

AP  (mm) 

2.5 

Water-fat  shift 

minimum 

Slice  thickness  (mm) 

2.5 

Shim 

auto 

Min.  TR/TE  (ms) 

1999  /  12 

Recon  voxel  size  (mm) 

2.5 

mDIXON 

no 

RF  avg  power  computed 
(W) 

1.371604 

Fold-over  suppression 

no 

Fat  suppression 

no 

SAR  /  head 

<  50  % 

Reconstruction  matrix 

96 

Water  suppression 

no 

Whole  body  /  level 

<  0.1  W/kg/ 
normal 

SENSE 

yes 

MTC 

no 

P  reduction  (AP) 

2.8 

Research  prepulse 

no 

B1  rms 

0.93  uT 

P  os  factor 

1 

Diffusion  mode 

no 

PNS  /  level  //  VUIIS  : 
dortch  : 

59  %  /  normal 

k-t  BLAST 

no 

Elastography  mode 

no 

Stacks 

1 

SAR  mode 

low 

Sound  Pressure  Level 
(dB) 

28.01657 

type 

parallel 

B1  mode 

default 

slices 

46 

SAR  Patient  data 

auto 

MOTION 

slice  gap 

user  defined 

PNS  mode 

low 

Cardiac  synchronization 

no 

gap  (mm) 

0 

Gradient  mode 

full  control 

Heart  rate  >  250  bpm 

no 

slice  orientation 

transverse 

max  strength 
(mT/m) 

33 

Respiratory 

compensation 

no 

fold -over  direction 

AP 

fat  shift  direction 

P 

max  slew  rate 
(T/m/s) 

130 

Navigator  respiratory 
comp 

no 

Minimum  number  of 
packages 

1 

Flow  compensation 

no 

Slice  scan  order 

default 

Temporal  slice  spacing 

equidistant 

PlanAlign 

no 

fMRI  echo  stabilisation 

no 

REST  slabs 

0 

NSA 

1 

Catheter  tracking 

no 

DYN/ANG 

Interactive  positioning 

no 

Angio  /  Contrast  enh. 

no 

Allow  table  movement 

no 

Quantitative  flow 

no 

nccr/  am 

i  n 

Manual  start 

yes 

Stacks 

i 

+Abuse  dynamic  loop 

no 

Stack  Offc.  AP 
(P=+mm) 

-5.532147 

Dynamic  study 

individual 

dyn  scans 

250 

RL  (L=+mm) 

-3.434925 

recon  multiplier 

1 

FH  (H  =+mm) 

35.45897 

dyn  scan  times 

shortest 

Ang.  AP  (deg) 

0 

FOV  time  mode 

default 

RL  (deg) 

-5.457352 

dummy  scans 

5 

FH  (deg) 

0 

immediate 

subtraction 

no 

fast  next  scan 

no 

synch,  ext.  device 

yes 

start  at  dyn. 

1 

interval  (dyn) 

119 

dyn  stabilization 

no 

prospect,  motion 

corr. 

no 

Keyhole 

no 

Arterial  Spin  labeling 

no 

POST/  PROC 

Preparation  phases 

full 

Interactive  FO 

no 

SENSE  ref.  scan 

no 

SmartPlan  survey 

no 

BO  field  map 

no 

B1  field  map 

no 

MIP/MPR 

no 

Images 

M,  no,  no,  no 

Autoview  image 

M 

Calculated  images 

no,  no,  no,  no 

Reference  tissue 

Grey  matter 

EPI  2D  phase  correction 

no 

Preset  window  contrast 

soft 

Reconstruction  mode 

real  time 

reuse  memory 

no 

Save  raw  data 

no 

Hardcopy  protocol 

no 

Ringing  filtering 

default 

Geometry  correction 

default 
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Q  Hospital  (2)  |  O20141021  CESTfMRI  (13)  52:37.4  |  E  FMRI_TRAI  LMAKI NG  04:14.0 


1 NFO  PAGE 

GEOMETRY 

CONTRAST 

Total  scan  duration 

04:14.0 

Multi-transmit 

no 

Scan  type 

Imaging 

Rel.  signal  level  (%) 

100 

Nucleus 

HI 

Scan  mode 

MS 

Act.  TR/TE  (ms) 

2000  /  25 

Coil  selection  1 

RX-lntf-1 

technique 

FFE 

Dyn.  scan  time 

00:02.000 

Xmit  Coil  selection 

MTX-Volume- 

T/R 

+  ZOOM 

no 

Time  to  kO 

00:15.0 

Contrast  enhancement 

no 

ACQ  matrix  M  x  P 

96x95 

User  def  elem  sel 

no 

Acquisition  mode 

cartesian 

ACQ  voxel  MPS  (mm) 

2.50/  2.50/ 

2.50 

element  selection 

All 

Fast  Imaging  mode 

EPI 

connection 

conn-A 

shot  mode 

single-shot 

REC  voxel  MPS  (mm) 

2.50  /  2.50  / 

2.50 

Coil  selection  2 

RX-lntf-2 

Echoes 

1 

element  selection 

All 

partial  echo 

no 

Scan  percentage  (%) 

100 

Dual  coil 

yes 

shifted  echo 

no 

Packages 

1 

CLEAR 

yes 

TE 

user  defined 

Min.  slice  gap  (mm) 

0 

body  tuned 

yes 

(ms) 

25 

EPI  factor 

37 

FOV  RL  (mm) 

240 

Flip  angle  (deg) 

63 

Act.  WFS  (pix)  /  BW 
(Hz) 

21.825  /  46.4 

AP  (mm) 

240 

TR 

user  defined 

FH  (mm) 

115 

(ms) 

2000 

BW  in  EPI  freq.  dir.  (Hz) 

2878.9 

Voxel  size  RL  (mm) 

2.5 

Ha  If  scan 

no 

Min.  WFS  (pix)  /  Max. 

BW  (Hz) 

21.770  /  46.5 

AP  (mm) 

2.5 

Water-fat  shift 

minimum 

Slice  thickness  (mm) 

2.5 

Shim 

auto 

Min.  TR/TE  (ms) 

1999  /  12 

Recon  voxel  size  (mm) 

2.5 

mDIXON 

no 

RF  avg  power  computed 
(W) 

1.371604 

Fold-over  suppression 

no 

Fat  suppression 

no 

SAR  /  head 

<  50  % 

Reconstruction  matrix 

96 

Water  suppression 

no 

Whole  body  /  level 

<  0.1  W/kg/ 
normal 

SENSE 

yes 

MTC 

no 

P  reduction  (AP) 

2.8 

Research  prepulse 

no 

B1  rms 

0.93  uT 

P  os  factor 

1 

Diffusion  mode 

no 

PNS  /  level  //  VUIIS  : 
dortch  : 

59  %  /  normal 

k-t  BLAST 

no 

Elastography  mode 

no 

Stacks 

1 

SAR  mode 

low 

Sound  Pressure  Level 
(dB) 

28.01657 

type 

parallel 

B1  mode 

default 

slices 

46 

SAR  Patient  data 

auto 

MOTION 

slice  gap 

user  defined 

PNS  mode 

low 

Cardiac  synchronization 

no 

gap  (mm) 

0 

Gradient  mode 

full  control 

Heart  rate  >  250  bpm 

no 

slice  orientation 

transverse 

max  strength 
(mT/m) 

33 

Respiratory 

compensation 

no 

fold -over  direction 

AP 

fat  shift  direction 

P 

max  slew  rate 
(T/m/s) 

130 

Navigator  respiratory 
comp 

no 

Minimum  number  of 
packages 

1 

Flow  compensation 

no 

Slice  scan  order 

default 

Temporal  slice  spacing 

equidistant 

PlanAlign 

no 

fMRI  echo  stabilisation 

no 

REST  slabs 

0 

NSA 

1 

Catheter  tracking 

no 

DYN/ANG 

Interactive  positioning 

no 

Angio  /  Contrast  enh. 

no 

Allow  table  movement 

no 

Quantitative  flow 

no 

nccr/  am 

i  n 

Manual  start 

yes 

Stacks 

i 

+Abuse  dynamic  loop 

no 

Stack  Offc.  AP 
(P=+mm) 

-5.532147 

Dynamic  study 

individual 

dyn  scans 

120 

RL  (L=+mm) 

-3.434925 

recon  multiplier 

1 

FH  (H  =+mm) 

35.45897 

dyn  scan  times 

shortest 

Ang.  AP  (deg) 

0 

FOV  time  mode 

default 

RL  (deg) 

-5.457352 

dummy  scans 

5 

FH  (deg) 

0 

immediate 

subtraction 

no 

fast  next  scan 

no 

synch,  ext.  device 

yes 

start  at  dyn. 

1 

interval  (dyn) 

119 

dyn  stabilization 

no 

prospect,  motion 

corr. 

no 

Keyhole 

no 

Arterial  Spin  labeling 

no 

POST/  PROC 

Preparation  phases 

full 

Interactive  FO 

no 

SENSE  ref.  scan 

no 

SmartPlan  survey 

no 

BO  field  map 

no 

B1  field  map 

no 

MIP/MPR 

no 

Images 

M,  no,  no,  no 

Autoview  image 

M 

Calculated  images 

no,  no,  no,  no 

Reference  tissue 

Grey  matter 

EPI  2D  phase  correction 

no 

Preset  window  contrast 

soft 

Reconstruction  mode 

real  time 

reuse  memory 

no 

Save  raw  data 

no 

Hardcopy  protocol 

no 

Ringing  filtering 

default 

Geometry  correction 

default 
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Q  Hospital  (2)  |  O20141021  CESTfMRI  (13)  52:37.4  |  E  FMRI_nback  08:30.0 


1 NFO  PAGE 

GEOMETRY 

CONTRAST 

Total  scan  duration 

08:30.0 

Multi-transmit 

no 

Scan  type 

Imaging 

Rel.  signal  level  (%) 

100 

Nucleus 

HI 

Scan  mode 

MS 

Act.  TR/TE  (ms) 

2500  /  25 

Coil  selection  1 

RX-lntf-1 

technique 

FFE 

Dyn.  scan  time 

00:02.500 

Xmit  Coil  selection 

MTX-Volume- 

T/R 

+  ZOOM 

no 

Time  to  kO 

00:18.7 

Contrast  enhancement 

no 

ACQ  matrix  M  x  P 

96x95 

User  def  elem  sel 

no 

Acquisition  mode 

cartesian 

ACQ  voxel  MPS  (mm) 

2.50/  2.50/ 

2.50 

element  selection 

All 

Fast  Imaging  mode 

EPI 

connection 

conn-A 

shot  mode 

single-shot 

REC  voxel  MPS  (mm) 

2.50  /  2.50  / 

2.50 

Coil  selection  2 

RX-lntf-2 

Echoes 

1 

element  selection 

All 

partial  echo 

no 

Scan  percentage  (%) 

100 

Dual  coil 

yes 

shifted  echo 

no 

Packages 

1 

CLEAR 

yes 

TE 

user  defined 

Min.  slice  gap  (mm) 

0 

body  tuned 

yes 

(ms) 

25 

EPI  factor 

37 

FOV  RL  (mm) 

240 

Flip  angle  (deg) 

63 

Act.  WFS  (pix)  /  BW 
(Hz) 

21.825  /  46.4 

AP  (mm) 

240 

TR 

user  defined 

FH  (mm) 

115 

(ms) 

2500 

BW  in  EPI  freq.  dir.  (Hz) 

2878.9 

Voxel  size  RL  (mm) 

2.5 

Ha  If  scan 

no 

Min.  WFS  (pix)  /  Max. 

BW  (Hz) 

21.770  /  46.5 

AP  (mm) 

2.5 

Water-fat  shift 

minimum 

Slice  thickness  (mm) 

2.5 

Shim 

auto 

Min.  TR/TE  (ms) 

1999  /  12 

Recon  voxel  size  (mm) 

2.5 

mDIXON 

no 

RF  avg  power  computed 
(W) 

1.097284 

Fold-over  suppression 

no 

Fat  suppression 

no 

SAR  /  head 

<  40  % 

Reconstruction  matrix 

96 

Water  suppression 

no 

Whole  body  /  level 

<  0.1  W/kg/ 
normal 

SENSE 

yes 

MTC 

no 

P  reduction  (AP) 

2.8 

Research  prepulse 

no 

B1  rms 

0.84  uT 

P  os  factor 

1 

Diffusion  mode 

no 

PNS  /  level  //  VUIIS  : 
dortch  : 

59  %  /  normal 

k-t  BLAST 

no 

Elastography  mode 

no 

Stacks 

1 

SAR  mode 

low 

Sound  Pressure  Level 
(dB) 

27.13099 

type 

parallel 

B1  mode 

default 

slices 

46 

SAR  Patient  data 

auto 

MOTION 

slice  gap 

user  defined 

PNS  mode 

low 

Cardiac  synchronization 

no 

gap  (mm) 

0 

Gradient  mode 

full  control 

Heart  rate  >  250  bpm 

no 

slice  orientation 

transverse 

max  strength 
(mT/m) 

33 

Respiratory 

compensation 

no 

fold -over  direction 

AP 

fat  shift  direction 

P 

max  slew  rate 
(T/m/s) 

130 

Navigator  respiratory 
comp 

no 

Minimum  number  of 
packages 

1 

Flow  compensation 

no 

Slice  scan  order 

default 

Temporal  slice  spacing 

equidistant 

PlanAlign 

no 

fMRI  echo  stabilisation 

no 

REST  slabs 

0 

NSA 

1 

Catheter  tracking 

no 

DYN/ANG 

Interactive  positioning 

no 

Angio  /  Contrast  enh. 

no 

Allow  table  movement 

no 

Quantitative  flow 

no 

nccr/  am 

i  n 

Manual  start 

yes 

Stacks 

i 

+Abuse  dynamic  loop 

no 

Stack  Offc.  AP 
(P=+mm) 

-5.532147 

Dynamic  study 

individual 

dyn  scans 

197 

RL  (L=+mm) 

-3.434925 

recon  multiplier 

1 

FH  (H  =+mm) 

35.45897 

dyn  scan  times 

shortest 

Ang.  AP  (deg) 

0 

FOV  time  mode 

default 

RL  (deg) 

-5.457352 

dummy  scans 

5 

FH  (deg) 

0 

immediate 

subtraction 

no 

fast  next  scan 

no 

synch,  ext.  device 

yes 

start  at  dyn. 

1 

interval  (dyn) 

119 

dyn  stabilization 

no 

prospect,  motion 

corr. 

no 

Keyhole 

no 

Arterial  Spin  labeling 

no 

POST/  PROC 

Preparation  phases 

full 

Interactive  FO 

no 

SENSE  ref.  scan 

no 

SmartPlan  survey 

no 

BO  field  map 

no 

B1  field  map 

no 

MIP/MPR 

no 

Images 

M,  no,  no,  no 

Autoview  image 

M 

Calculated  images 

no,  no,  no,  no 

Reference  tissue 

Grey  matter 

EPI  2D  phase  correction 

no 

Preset  window  contrast 

soft 

Reconstruction  mode 

real  time 

reuse  memory 

no 

Save  raw  data 

no 

Hardcopy  protocol 

no 

Ringing  filtering 

default 

Geometry  correction 

default 
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CD  Hospital  (2)  |  O20141021  CEST  fMRI  (13)  52:37.4  |  En_3D_TFE_quantGeo  00:55.7 


1 NFO  PAGE 

GEOMETRY 

CONTRAST 

Total  scan  duration 

00:55.7 

Multi-transmit 

no 

Scan  type 

Imaging 

Rel.  signal  level  (%) 

100 

Nucleus 

HI 

Scan  mode 

3D 

Act.  TR/TE  (ms) 

2.8/  1.44 

Coil  selection  1 

RX-lntf-1 

technique 

FFE 

ACQ  matrix  M  x  P 

192  x  190 

Xmit  Coil  selection 

MTX-Volume- 

+  ZOOM 

no 

ACQ  voxel  MPS  (mm) 

1.25  /  1.26  / 

T/R 

Contrast  enhancement 

Tl 

2.50 

User  def  elem  sel 

no 

Acquisition  mode 

cartesian 

REC  voxel  MPS  (mm) 

0.94/  0.94/ 

element  selection 

All 

Fast  Imaging  mode 

TFE 

2.50 

connection 

conn-A 

3D  non -selective 

no 

Scan  percentage  (%) 

98.9899 

Coil  selection  2 

RX-lntf-2 

shot  mode 

multishot 

TFE  shots 

13 

element  selection 

All 

TFE  factor 

256 

TFE  dur.  shot  /  acq  (ms) 

1675.7  /  729.0 

Dual  coil 

yes 

3D  free  factor 

no 

Min.  Tl  delay 

386.7124 

CLEAR 

yes 

startup  echoes 

default 

Act.  WFS  (pix)  /  BW 

0.778  /  1302.1 

body  tuned 

yes 

+TFE  followup  echoes 

0 

(Hz) 

FOV  RL  (mm) 

240 

shot  interval 

user  defined 

Min.  WFS  (pix)  /  Max. 

BW  (Hz) 

0.774  /  1308.4 

AP  (mm) 

240 

(ms) 

4500 

RF  avg  power  computed 
(W) 

0.8993402 

FH  (mm) 

165 

profile  order 

linear 

Voxel  size  RL  (mm) 

1.25 

turbo  direction 

radial 

SAR  /  head 

<  33  % 

AP  (mm) 

1.254902 

CENTRA  (spiral) 

no 

Whole  body  /  level 

0.0  W/kg  / 

FH  (mm) 

2.5 

Echoes 

1 

normal 

Recon  voxel  size  (mm) 

0.94 

partial  echo 

no 

B1  rms 

0.76  uT 

Fold-over  suppression 

no 

shifted  echo 

no 

PNS  /  level  //  VUIIS  : 

59  %  /  normal 

Slice  oversampling 

default 

TE 

shortest 

dortch  : 

RF  select.  FOS 

no 

Flip  angle  (deg) 

7 

Sound  Pressure  Level 
(dB) 

29.78219 

Reconstruction  matrix 

256 

TR 

shortest 

SENSE 

yes 

Ha  If  scan 

no 

MOTI  or 

\l 

P  reduction  (AP) 

2 

Water -fat  shift 

minimum 

Cardiac  synchronization 

no 

P  os  factor 

1 

Shim 

auto 

Heart  rate  >  250  bpm 

no 

S  reduction  (FH) 

2 

mDIXON 

no 

Respiratory 

compensation 

no 

k-t  BLAST 

no 

Fat  suppression 

no 

Overcontiguous  slices 

no 

Water  suppression 

no 

Navigator  respiratory 
comp 

Stacks 

1 

TFE  prepulse 

invert 

Flow  compensation 

no 

slices 

66 

slice  selection 

no 

fMRI  echo  stabilisation 

no 

slice  orientation 

transverse 

delay 

user  defined 

Motion  smoothing 

no 

fold -over  direction 

AP 

(ms) 

1300 

NSA 

1 

fat  shift  direction 

L 

PSIR 

no 

DYN/  ANi 

G 

Chunks 

1 

+inv  pulse  type 

+B1  opt  (low 
BW) 

PlanAlign 

no 

Angio  /  Contrast  enh. 

no 

REST  slabs 

0 

MTC 

no 

Quantitative  flow 

no 

Catheter  tracking 

no 

T2prep 

no 

CENTRA 

no 

Interactive  positioning 

no 

Research  prepulse 

no 

Manual  start 

no 

Allow  table  movement 

no 

Diffusion  mode 

no 

+Abuse  dynamic  loop 

no 

Elastography  mode 

no 

Dynamic  study 

no 

urrw  anu 

Stacks 

SAR  mode 

high 

Arterial  Spin  labeling 

no 

1 

Stack  Offc.  AP 

-4.247866 

POST/  PROC 

di  moae 

aeTdUit 

(P=+mm) 

SAR  Patient  data 

auto 

Preparation  phases 

auto 

RL  (L=+mm) 

-3.434925 

PNS  mode 

low 

no 

i  nteractive  fu 

FH  (H  =+mm) 

35.88707 

Gradient  mode 

full  control 

CCMCC  rr\f 

no 

btNbt  rGT.  scan 

Ang.  AP  (deg) 

o 

max  strength 

33 

SmartPlan  survey 

no 

RL  (deg) 

-5.457352 

(mT/m) 

BO  field  map 

no 

FH  (deg) 

0 

max  slew  rate 

166 

B1  field  map 

no 

(T/m/s) 

MIP/MPR 

no 

Images 

M,  no,  no,  no 

Autoview  image 

M 

Calculated  images 

no,  no,  no,  no 

Reference  tissue 

Grey  matter 

Preset  window  contrast 

soft 

Reconstruction  mode 

immediate 

Save  raw  data 

no 

Hardcopy  protocol 

no 

Ringing  filtering 

rectangular 

Geometry  correction 

default 

Elliptical  k-space  shutter 

default 
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1 NFO  PAGE 

GEOMETRY 

CONTRAST 

Total  scan  duration 

01:02.7 

Multi-transmit 

no 

Scan  type 

Imaging 

Rel.  signal  level  (%) 

100 

Nucleus 

HI 

Scan  mode 

3D 

Act.  TR/TEl/delta  TE 
(ms) 

34/  3.3/  3.2 

Coil  selection  1 

RX-lntf-1 

technique 

FFE 

Xmit  Coil  selection 

MTX-Volume- 

T/R 

loop  order 

zy  order 

ACQ  matrix  M  x  P 

240  x  240 

+  ZOOM 

no 

ACQ  voxel  MPS  (mm) 

1.00/  1.00/ 

5.00 

User  def  elem  sel 

no 

Contrast  enhancement 

T1 

element  selection 

All 

Acquisition  mode 

cartesian 

REC  voxel  MPS  (mm) 

0.94/  0.94/ 

5.00 

connection 

conn-A 

Fast  Imaging  mode 

none 

Coil  selection  2 

RX-lntf-2 

3D  non-selective 

no 

Scan  percentage  (%) 

100 

element  selection 

All 

Echoes 

10 

Act.  WFS  (pix)  /  BW 
(Hz) 

1.418/  714.8 

Dual  coil 

yes 

partial  echo 

no 

Min.  WFS  (pix)  /  Max. 

BW  (Hz) 

1.125  /  900.9 

CLEAR 

yes 

shifted  echo 

no 

body  tuned 

yes 

TE  first 

shortest 

RF  avg  power  computed 
(W) 

0.1353424 

FOV  AP  (mm) 

240 

echospacing 

shortest 

RL  (mm) 

240 

flyback 

yes 

SAR  /  head 

<  5  % 

FH  (mm) 

60 

Flip  angle  (deg) 

8 

Whole  body  /  level 

0.0  W/kg  / 
normal 

Voxel  size  AP  (mm) 

1 

TR 

shortest 

RL  (mm) 

1 

Ha  If  scan 

no 

B1  rms 

0.29  uT 

FH  (mm) 

5 

Water -fat  shift 

user  defined 

PNS  /  level  //  VUIIS  : 
dortch  : 

60  %  /  normal 

Recon  voxel  size  (mm) 

0.9375 

(pixels) 

1.4 

Fold-over  suppression 

no 

Shim 

PB-volume 

Sound  Pressure  Level 
(dB) 

29.7646 

Slice  oversampling 

default 

ShimAlign 

no 

RF  select.  FOS 

no 

mDIXON 

no 

MOTION 

Reconstruction  matrix 

256 

Fat  suppression 

no 

Cardiac  synchronization 

no 

SENSE 

yes 

Water  suppression 

no 

Heart  rate  >  250  bpm 

no 

P  reduction  (RL) 

2 

MTC 

no 

Respiratory 

compensation 

no 

P  os  factor 

1 

Research  prepulse 

no 

S  reduction  (FH) 

1 

Diffusion  mode 

no 

Navigator  respiratory 
comp 

no 

k-t  BLAST 

no 

Elastography  mode 

no 

Flow  compensation 

yes 

Overcontiguous  slices 

no 

SAR  mode 

low 

fMRI  echo  stabilisation 

no 

Stacks 

1 

B1  mode 

default 

NSA 

1 

slices 

12 

SAR  Patient  data 

auto 

DYN/AN< 

G 

slice  orientation 

transverse 

PNS  mode 

low 

Angio  /  Contrast  enh. 

inflow 

fold -over  direction 

RL 

Gradient  mode 

maximum 

Quantitative  flow 

no 

fat  shift  direction 

P 

SofTone  mode 

no 

Tone  pulse 

no 

Chunks 

1 

PlanAlign 

no 

Manual  start 

no 

+Abuse  dynamic  loop 

no 

REST  slabs 

0 

Dynamic  study 

no 

Catheter  tracking 

no 

Arterial  Spin  labeling 

no 

Interactive  positioning 

no 

Allow  table  movement 

no 

POST/  PROC 

Preparation  phases 

auto 

OFFC/AN 

IG 

Interactive  FO 

no 

Stacks 

1 

SENSE  ref.  scan 

no 

Stack  Offc.  AP 
(P=+mm) 

-5.532147 

SmartPlan  survey 

no 

RL  (L=+mm) 

-3.434925 

BO  field  map 

no 

FH  (H  =+mm) 

35.45897 

B1  field  map 

no 

Ang.  AP  (deg) 

0 

MIP/MPR 

no 

Images 

M,  R,  1,  no 

RL  (deg) 

-5.457352 

FH  (deg) 

0 

Autoview  image 

M 

Shim  Size  AP  (mm) 

180.7495 

Calculated  images 

T2,  no,  no,  no 

RL  (mm) 

134.1432 

T2*  clipvalue  (ms) 

100 

FH  (mm) 

114.2866 

Reference  tissue 

Grey  matter 

Offc.  AP  (P=+mm) 

-1.87074 

Preset  window  contrast 

soft 

Reconstruction  mode 

immediate 

RL  (L=+mm) 

-3.660397 

FH  (H  =+mm) 

33.86363 

Save  raw  data 

no 

Ang.  AP  (deg) 

2.26406 

Hardcopy  protocol 

no 

RL  (deg) 

-6.580403 

Ringing  filtering 

default 

FH  (deg) 

4.637685 

Geometry  correction 

default 

Elliptical  k-space  shutter 

no 
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Quantitative  magnetization  transfer  (qMT)  imaging  yields  indices  describing  the  interactions  between  free  water 
protons  and  immobile  macromolecular  protons.  These  indices  include  the  macromolecular  to  free  pool  size  ratio 
(PSR),  which  has  been  shown  to  be  correlated  with  myelin  content  in  white  matter.  Because  of  the  long  scan 
times  required  for  whole-brain  imaging  (=  20-30  min),  qMT  studies  of  the  human  brain  have  not  found  wide¬ 
spread  application.  Herein,  we  investigated  whether  the  increased  signal-to-noise  ratio  available  at  7.0  T  could 
be  used  to  reduce  qMT  scan  times.  More  specifically,  we  developed  a  selective  inversion  recovery  (SIR)  qMT  imag¬ 
ing  protocol  with  a  i)  novel  transmit  radiofrequency  (B^  )  and  static  field  (B0)  insensitive  inversion  pulse,  ii)  turbo 
field-echo  readout,  and  iii)  reduced  TR.  In  vivo  qMT  data  were  obtained  in  the  brains  of  healthy  volunteers  at  7.0  T 
using  the  resulting  protocol  (scan  time=  40  s/slice,  resolution  =2  x  2^3  mm3).  Reliability  was  also  assessed  in  re¬ 
peated  acquisitions.  The  results  of  this  study  demonstrate  that  SIR  qMT  imaging  can  be  reliably  performed  within 
the  radiofrequency  power  restrictions  present  at  7.0  T,  even  in  the  presence  of  large  and  &  inhomogeneities. 
Consistent  with  qMT  studies  at  lower  field  strengths,  the  observed  PSR  values  were  higher  in  white  matter 
(mean  ±  SD  =  17.6  ±  1.3%)  relative  to  gray  matter  (10.3  ±  1.6%)  at  7.0  T.  In  addition,  regional  variations  in  PSR 
were  observed  in  white  matter.  Together,  these  results  suggest  that  qMT  measurements  are  feasible  at  7.0  T  and 
may  eventually  allow  for  the  high-resolution  assessment  of  changes  in  composition  throughout  the  normal  and 
diseased  human  brain  in  vivo. 

©  2012  Elsevier  Inc.  All  rights  reserved. 


Introduction 

In  addition  to  the  free  water  protons  typically  observed  in  mag¬ 
netic  resonance  imaging  (MRI),  there  are  protons  residing  on  immo¬ 
bile  macromolecules  in  tissue  (Wolff  and  Balaban,  1989).  Typical 
imaging  sequences  do  not  directly  detect  this  pool  of  protons  because 
they  exhibit  very  short  transverse  relaxation  times  (~  10  !Js)  and,  there¬ 
fore,  lose  coherence  before  their  signal  can  be  captured.  This  macromol¬ 
ecule  proton  pool  can,  however,  be  indirectly  detected  by  exploiting  its 
interactions  with  the  free  water  pool  via  chemical  exchange  and/or  dipo¬ 
lar  mechanisms  [referred  to  together  as  the  magnetization  transfer  (MT) 
effect].  Previous  phantom  studies  (Koenig,  1991;  Kucharczyk  et  al.,  1994) 
have  shown  that  the  bulk  of  the  MT  effect  in  white  matter  (WM)  arises 
from  myelin-associated  lipids,  which  suggests  that  MT  contrast  may  be 
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a  more  specific  marker  for  myelin  pathology  than  conventional  imaging 
methods.  As  a  result,  there  is  considerable  interest  in  exploiting  MT  con¬ 
trast  to  assay  changes  in  myelination  associated  with  a  number  of  dis¬ 
eases  [e.g.,  multiple  sclerosis  (Catalaa  et  al.,  2000;  Filippi  and  Rocca, 
2004;  Gass  et  al.,  1994;  Kalkers  et  al.,  2001)  and  neuropsychiatric  dis¬ 
eases  (Bruno  et  al.,  2004;  Kabani  et  al.,  2002a,  2002b)]. 

MT  contrast  can  be  generated  by  applying  an  off-resonance 
radiofrequency  (RF)  prepulse  to  selectively  saturate  the  spectrally 
broad  macromolecular  proton  pool  (Wolff  and  Balaban,  1989).  This 
saturation  then  transfers  to  the  free  water  proton  pool  via  MT, 
resulting  in  a  decrease  in  the  observed  free  water  signal.  The  magni¬ 
tude  of  this  effect  can  be  characterized  by  a  semi-quantitative  metric 
known  as  the  magnetization  transfer  ratio  (Dousset  et  al.,  1992): 
MTR=  1  -  Ssat/S0,  where  Ssat  and  S0  are  the  observed  signal  intensi¬ 
ties  with  and  without  the  application  of  an  MT  saturation  prepulse, 
respectively.  Although  the  MTR  has  been  shown  to  correlate  with  my¬ 
elin  content  (Odrobina  et  al.,  2005;  Schmierer  et  al.,  2004),  it  is  also  sen¬ 
sitive  to  the  choice  of  experimental  parameters  such  as  RF  power  (Berry 
et  al.,  1999)  as  well  as  non-MT-specific  NMR  parameters  such  as  tis¬ 
sue  relaxation  times  (Henkelman  et  al.,  1993).  As  a  result,  quantitative 
MT  (qMT)  approaches  have  been  developed.  These  qMT  approaches 
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quantify  distinct  tissue  characteristics  (e.g.,  the  size  of  the  macromolecu- 
lar  pool,  rate  of  MT  exchange)  rather  than  the  combined  effect  of  multi¬ 
ple  tissue  and/or  acquisition  parameters.  As  such,  qMT  measures  are 
thought  to  yield  more  specific  information  on  tissue  composition  than 
the  MTR. 

Pulsed  saturation  qMT  imaging  (Graham  and  Henkelman,  1997; 
Pike,  1996;  Sled  and  Pike,  2000,  2001)  has  received  considerable  at¬ 
tention  for  application  in  humans  in  vivo  because  it  allows  for  the 
rapid  collection  of  qMT  data  within  the  hardware  constraints  of 
most  clinical  systems.  This  approach  involves  a  steady-state,  spoiled 
gradient-echo  acquisition  interleaved  with  an  MT-preparation  pulse. 
By  collecting  images  over  a  range  of  MT  pulse  offset  frequencies 
and/or  powers  and  fitting  the  resulting  data  to  a  two-pool  model  of 
the  MT  effect,  one  can  extract  parameters  such  as  the  macromolecular 
to  free  pool  size  ratio  (PSR)  and  the  rate  of  MT  exchange.  Previous 
work  has  shown  that  the  PSR  is  correlated  with  myelin  content 
(Odrobina  et  al.,  2005;  Ou  et  al.,  2009;  Schmierer  et  al.,  2007; 
Underhill  et  al.,  2011).  The  relationship  between  the  rate  of  MT  ex¬ 
change  and  underlying  tissue  composition  is  less  clear;  however, 
previous  work  has  suggested  that  the  rate  of  MT  exchange  may  re¬ 
flect  changes  within  the  myelin  lipid  structure  (Smith  et  al.,  2009). 

Unfortunately,  qMT  imaging  has  not  found  widespread  applica¬ 
tion  in  practice.  This  can  be  attributed  in  part  to  the  long  scan  times 
(~  20-30  min  for  whole-brain  imaging)  required  to  collect  images 
at  multiple  offset  frequencies  and/or  powers.  The  number  of  total  im¬ 
ages  required  can  be  reduced  by  designing  optimal  sampling  strate¬ 
gies  (Cercignani  and  Alexander,  2006;  Levesque  et  al.,  2011)  or  by 
fixing  certain  model  parameters  in  the  fitting  procedure  (Underhill 
et  al.,  2009,  2011).  Potentially  more  efficient  strategies  based  upon 
steady-state  free-precession  (SSFP)  sequences  (Garcia  et  al.,  2010; 
Gloor  et  al.,  2008)  may  also  be  employed. 

As  an  alternative,  or  perhaps  in  combination  with  these  strategies, 
one  could  translate  qMT  imaging  approaches  to  higher  field  strengths. 
The  resulting  increase  in  SNR  could  then  be  used  to  obtain  more  reliable 
estimates  of  MT  parameters  or  traded  to  reduce  scan  times  and/or  in¬ 
crease  resolution.  To  date,  qMT  studies  in  humans  in  vivo  have  been  pri¬ 
marily  limited  to  1.5  and  3.0  T  and  we  are  aware  of  only  one  report 
(Mougin  et  al.,  2010)  of  MT  parameters  in  the  human  brain  in  vivo  at 
7.0  T.  The  translation  of  pulsed  saturation  approaches  to  7.0  T  faces 
two  primary  challenges:  i)  RF  power  limitations  [e.g.,  specific  ab¬ 
sorption  ratio  (SAR)  limitations]  and  ii)  transmit  RF  (^_H)  and  static 
magnetic  field  (B0)  inhomogeneities.  SSFP-based  approaches  may 
also  be  limited  at  high  field  by  banding  artifacts  associated  with  B0 
inhomogeneities.  In  contrast,  selective  inversion  recovery  (SIR)  qMT 

imaging  (Edzes  and  Samulski,  1977;  Gochberg  et  al.,  1997),  which  is 
based  upon  measuring  the  biexponential  recovery  of  the  free  water 
pool  in  the  presence  of  MT  after  an  on-resonance  inversion  pulse,  has 
been  suggested  (Dortch  et  al.,  201 1)  to  be  less  sensitive  to  these  issues. 
Note  that  this  approach  is  similar  to  the  stimulated  echo  approach  pro¬ 
posed  by  Ropele  et  al.  (2003);  therefore,  both  approaches  may  be  well 
suited  for  qMT  imaging  at  7.0  T. 

In  this  study,  we  have  investigated  the  feasibility  of  using  the  SIR 
approach  for  high  field  qMT  imaging  of  the  human  brain.  More  specif¬ 
ically,  we  have  translated  our  previously  published  3.0-T  SIR  protocol 
(Dortch  et  al.,  201 1 )  to  7.0  T  with  two  significant  modifications.  First, 
we  incorporated  anovel  B^-and  !1  ^-insensitive  composite  inversion 
pulse  to  ensure  a  more  uniform  inversion  of  the  free  water  pool  over 
the  whole  brain.  Second,  we  transitioned  from  a  turbo-spin  echo 
readout  (TSE)  to  a  turbo  field-echo  readout  (TFE) — similar  to  an 
MP-RAGE  sequence  (Mugler  and  Brookeman,  1990)— as  the  former 
is  susceptible  to  Bf -related  artifacts  (due  to  imperfect  refocusing) 
and  is  SAR-limited  at  high  field.  The  TFE  readout  has  the  added  bene¬ 
fit  of  covering  k-space  more  efficiently  than  the  TSE  readout,  which,  in 
combination  with  some  additional  protocol  optimization,  allowed  us 
to  transition  from  a  single-slice  approach  at  3.0  T  to  a  whole-brain  ap¬ 
proach  at  7.0  T  (~  40  s/slice  at  2.0 *  2.0 x  3.0  mm3  resolution).  Using 


this  protocol,  in  vivo  qMT  data  were  obtained  in  the  brains  of  13 
healthy  volunteers  at  7.0  T.  To  assess  the  reproducibility  of  the  tech¬ 
nique,  six  of  the  healthy  volunteers  were  scanned  twice.  Additional 
numerical  simulations  were  performed  to  determine  the  effect  of 
TFE  readout  on  our  qMT  parameter  maps. 


Theory 


Consider  free  water  (f)  and  macromolecular  (m)  proton  pools  be¬ 
tween  which  MT  can  occur.  Define  unique  equilibrium  magnetiza¬ 
tions  (M0f  and  M0m),  spin-lattice  relaxation  rates  (R1f  and  R1m),  and 
spin-spin  relaxation  rates  (R2f  and  R2m)  for  each  pool  as  well  as  an 
MT  rate  from  the  macromolecular  to  the  free  pool  (kmf)— the  rate  in 
the  other  direction  can  be  determined  from  kfm  =  kmfM0m/ M0f.  As¬ 
sume  MT  of  transverse  magnetization  to  be  negligible  because  of 
the  short  T2  of  the  macromolecular  pool.  In  this  case,  the  transverse 
components  of  the  macromolecular  pool  can  be  ignored.  The  time 
evolution  of  the  remaining  x,  y,  and  z  components  of  the  magnetiza¬ 
tion  vector  M  =  [Mxf  Myf  M^  Mzm]T  during  a  constant  amplitude  RF 
pulse  can  be  expressed  in  matrix  form  as  (Portnoy  and  Stanisz,  2007) 
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!1  (jo  is  the  frequency  offset  from  resonance  for  the  RF  pulse,  co-,  is  the 
frequency  of  precession  about  the  RF  pulse,  and  (J)  is  the  phase  of 
the  RF  pulse  in  the  transverse  plane.  The  standard  Bloch  equations 
implicitly  assume  a  Lorentzian  lineshape,  which  is  invalid  for  the 
macromolecular  proton  pool.  As  a  result,  the  Bloch  equations  for  the 
macromolecular  pool  have  been  replaced  in  Eq.  (2)  by  a  single  longi¬ 
tudinal  component  whose  saturation  is  governed  by  the  rate  Rrf  = 

TT^I  m  m 

2g  (  ou),  where  g  is  the  lineshape  function  of  the  macromolecu¬ 
lar  pool.  When  applying  off-resonance  irradiation,  a  super-Lorentzian 
lineshape  is  typically  used  to  model  biological  macromolecular  pro¬ 
tons  (Morrison  et  al.,  1995).  Because  the  super-Lorentzian  exhibits 
an  on-resonance  singularity,  Gaussian  (Gochberg  and  Gore,  2007)  or 
super-Lorentzian  functions  extrapolated  from  a  1  kHz  offset  (Gloor 
et  al.,  2008)  are  typically  used  to  model  the  macromolecular  pool 
lineshape  pool  during  on-resonance  irradiation. 

The  general  solution  to  this  system  of  equations  can  be  expressed 
as 


Mdtb  %  exp6AtbM60b  \d  '/zexpSAtb—  l]A  63b 

where  M(0)  is  the  initial  condition  of  the  system  and  I  is  an  identity 
matrix.  The  same  expression  can  be  used  to  describe  the  system  dur¬ 
ing  free  precession  (i.e. ,  when  cu-,  =  0).  In  this  case,  the  solution  can  be 
further  simplified  by  noting  that  the  z-component  is  decoupled  from 
the  x-  and  y-components,  resulting  in  the  following  expression  for  the 
longitudinal  magnetization  vector  Mz  =  [Mzf  Mzm]T 
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where  M0  =  [M0f  M0m]T  and  Az  is  the  lower-right  quadrant  of  A  with 
Rrf  =  0.  Expanding  the  matrix  exponentials  in  this  expression  yields 
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where  A-1-7-  are  the  negative  eigenvalues  of  Az  and  U  is  a  matrix 
whose  columns  are  the  corresponding  eigenvectors.  From  Eq.  (5),  it 
can  be  seen  that  recovers  as  a  biexponential  function  governed 
by  the  fast  and  slow  rate  constants  A-1-  and  A- ,  respectively,  during 
free  precession.  As  described  below,  one  can  obtain  estimates  of 
qMT  parameters  (e.g.,  PSR  and  kmf)  by  measuring  this  biexponential 
recovery. 

Methods 

Pulse  sequence 


lineshape),  T2f  =  60  ms,  kmf  =  15  s-  \  and  PSR  =  15%.  The  results  from 
these  simulations  indicate  that  this  pulse  train  saturates  both  pools 
[Mzf(tsat)/M0f<  0.01  and  Mzm(tsat)/M0m<  0.06]  over  the  range  of 
expected  B-1-  values  (B-1-  /B-1-  =  0.3-1 .0)  in  the  human  brain 

1  1  .actual  1, nominal 

at  7.0  T  [n.b.,  the  manufacturer-provided  power  optimization  tended 
to  yield  a  mean  K'actuJKwnMnaibl.O  (Moore  et  al.,  2010)]. 

Plugging  the  initial  condition  of  Mz(td  =  0)  =0  into  Eq.  (4),  signal 
equations  can  be  generated  for  the  predelay  period  of  the  SIR-TFE  se¬ 
quence.  The  ending  values  for  this  period  can  then  be  used  as  the  ini¬ 
tial  condition  for  the  inversion  recovery  period,  taking  account  for  the 
effect  of  the  inversion  pulse 

Mz3t|3t5  %  SMz6t-l=;  361= 

where  Sis  a  diagonal  matrix  with  elements  that  account  for  the  inver¬ 
sion  of  the  free  pool  (Sf=-  1  denotes  complete  inversion)  and  the 
saturation  of  the  macromolecular  pool  (Sm  =  1  denotes  no  saturation) 
and  t-1-7-  is  the  time  immediately  before/after  the  pulse.  This  yields 
the  final  expression  for  the  evolution  of  Mz  during  the  SIR  period  of 
the  sequence 


The  SIR  qMT  sequence  (Fig.  1 )  used  herein  is  similar  to  the  inver¬ 
sion  recover  sequence  used  to  measure  T1  with  two  modifications. 
First,  short  inversion  times  (~  10  ms  or  less)  are  sampled  in  order 

to  capture  the  fast-recovering  A-1-  component  of  the  biexponential  re¬ 
covery.  Second,  a  T2-selective  inversion  pulse  is  applied.  This  is 
achieved  via  a  low  power  inversion  pulse  whose  duration  is  much 
longer  than  the  T2  of  the  macromolecular  pool  (T2m  =  10  !Js)  and 
much  shorter  than  the  T2  of  the  free  water  pool  (T2f=  10-100  ms). 
Ideally,  this  pulse  inverts  M^  with  minimal  saturation  of  Mzm.  In 
other  words,  this  pulse  maximizes  the  difference  between  the  pools 
and,  in  turn,  the  sensitivity  of  the  signal  to  MT.  This  is  followed  by  a 
variable  duration  inversion  recovery  period  to  sample  the  transient 
biexponential  recovery  of  Mzf  and  a  center-out  TFE  readout  (SIR-TFE) 
to  efficiently  sample  k-space.  For  inversion  recovery  acquisitions,  a 
predelay  time  td  =  5/A-  is  commonly  employed  to  ensure  full  recov¬ 
ery  of  M^.  Flowever,  if  one  can  assume  that  the  longitudinal  magneti¬ 
zation  of  both  pools  is  approximately  zero  at  the  end  of  the  readout, 
the  effect  of  a  shorter  predelay  period  can  be  accounted  for  in  the  sig¬ 
nal  model,  allowing  one  to  reduce  td  (and  scan  times)  without  biasing 
the  estimated  parameters.  This  assumption  has  been  previously 
shown  to  hold  true  for  a  TSE  readout  (Gochberg  and  Gore,  2007); 
however,  this  cannot  be  assumed  for  the  TFE  readout  employed  here¬ 
in.  As  a  result,  we  empirically  designed  a  train  of  RF  pulses  [number  of 
pulses  =  32,  a  =  135°,  pulse  spacing  =  20  ms,  pulse  train  duration 
(tsat)  =  620  ms]  to  saturate  both  pools  following  the  TFE  readout.  To 
assess  the  effect  of  this  pulse  train  on  the  longitudinal  magnetization 
of  both  pools,  numerical  simulations  were  performed  via  Eq.  (3)  and 
the  following  parameters:  R1m  =  R1f  =  0.8  s-  \  T2m  =  10  !Js  (Gaussian 
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Fig.  1.  SIR-TFE  pulse  sequence  diagram.  The  sequence  employs  i)  a  composite  inversion 
pulse  (Fig.  2)  designed  to  uniformly  invert  Mzf  over  a  range  of  expected  !1B0and  Bi"*" 
values  with  minimal  macromolecular  pool  saturation,  ii)  a  variable  duration  inversion 
recovery  (SIR)  period  to  sample  the  free  pool  recovery,  iii)  a  TFE  readout  to  efficiently 
cover  k-space,  iv)  a  pulse  train  to  saturate  (SAT)  the  free  and  macromolecular  pools 
(allows  tdb 5/A- ),  and  v)  a  predelay  (PD)  period  to  allow  for  partial  Mz  recovery. 
Legend:  tj  =  inversion  time,  td  =  predelay,  t  =  TFE  pulse-to-pulse  interval,  ACQ  = 
acquisition. 
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In  addition  to  these  pulse  sequence  modifications,  a  novel 
64-element  composite  inversion  pulse  was  designed  and  employed 
herein  to  ensure  a  uniform  inversion  of  Mzf  over  the  range  of 
and  !1B0  values  previously  measured  in  the  human  brain  at  7.0  T 
(Moore  et  al.,  2010).  The  optimization  procedure  (Moore  et  al.,  2010) 
tended  to  produce  high  power  pulses  with  suboptimal  T2-selectivity. 
As  a  result,  we  included  an  additional  RF  power  constraint  into  the 
procedure,  which  was  weighted  against  the  uniform  inversion  con¬ 
straint.  The  resulting  amplitudes  and  phases  of  the  subpulses  are 
shown  in  Fig.  2.  To  evaluate  the  pulse's  performance,  Sf  and  Sm  were 
estimated  from  Eq.  (6)  by  propagating  Eq.  (3)  through  each  of  the  64 
subpulses  [neglecting  T1  relaxation  and  MT  during  the  pulse  and 
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Fig.  2.  Composite  inversion  pulse  amplitudes  (a),  phases  (b),  predicted  free  water  inversion 
efficiency  Sf  (c),  and  predicted  macromolecular  saturation  fractions  Sm  (d).  Sf=—  1  denotes 
complete  inversion;  Sm=  1  denotes  no  saturation.  The  two  zero-amplitude  discontinuities 
in  the  RF  pulse  (a)  are  a  consequence  of  the  power  constraint  used  in  the  minimization  pro¬ 
cedure.  The  RF  phase  (b)  of  the  pulse  at  these  discontinuities  is  arbitrary;  therefore,  the 
phase  at  these  points  was  set  based  upon  linear  interpolation  of  the  neighboring  RF  phases 
for  display  purposes. 
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assuming  a  Gaussian  macro  molecular  pool  lineshape  with  T2m  =  10  !Js 
(Gochberg  and  Gore,  2007)].  From  this  procedure,  the  pulse  is  predicted 
to  yield  a  uniform  inversion  of  over  a  wide  range  of  and  !1B0 
values  without  complete  saturation  of  Mzm. 


Numerical  simulations 

The  TFE  readout  employed  herein  effectively  blurs  the  image  along 
the  phase-encoding  direction  according  to  its  readout  point-spread 
function  (PSF),  which  is  a  complex  function  of  the  sequence  timings 
and  the  NMR  parameters  of  the  tissue  (Constable  and  Gore,  1 992).  If 
the  readout  PSF  is  constant  as  a  function  of  tj,  then  its  effect  will  be 
to  simply  blur  the  final  MT  parameter  maps.  If,  however,  the  readout 
PSF  changes  as  a  function  of  tj,  each  image  will  be  blurred  to  a  differ¬ 
ent  degree,  potentially  biasing  the  final  parameter  maps. 

To  evaluate  this  effect,  the  SIR-TFE  signal  arising  from  a  one¬ 
dimensional  (ID)  test  object  was  numerically  simulated.  As  shown 
in  Fig.  3,  MT  parameters  were  defined  for  test  object  regions 
representing  white  matter  (WM),  gray  matter  (GM),  and  cerebro¬ 
spinal  fluid  (CSF).  For  each  region  and  tj,  the  signal  evolution  during 
each  RF  pulse  and  precession  period  of  the  TFE  readout  was  simulat¬ 
ed  from  Eq.  (3)  with  the  imaging  parameters  in  the  Data  acquisition 
section— using  Mz(tj)  from  Eq.  (7)  as  the  initial  condition  and  re¬ 
placing  each  time-varying  excitation  pulse  with  a  constant  ampli¬ 
tude  pulse  of  equivalent  flip  angle  or  root-mean-squared  power 
(Ramani  et  al.,  2002)  for  the  free  water  or  macromolecular  pool,  re¬ 
spectively.  Complete  spoiling  of  transverse  magnetization  was  as¬ 
sumed  prior  to  each  RF  pulse.  The  resulting  Mzf  immediately  after 
each  RF  pulse  was  taken  to  represent  the  signal  as  a  function  of 
echo  number.  The  signal  was  then  re-ordered  to  account  for  the 
k-space  trajectory  and  SENSE  acceleration  used,  and  the  resulting 
re-ordered  signal  was  taken  to  represent  a  k-space  filter.  To  apply 
the  k-space  filters  to  the  ID  test  object,  each  uniform  object  region 
was  Fourier  transformed  into  k-space,  multiplied  by  its  correspond¬ 
ing  k-space  filter,  and  inverse  Fourier  transformed  back  into  image 
space.  The  resulting  object  regions  were  then  summed  to  generate 


voxel  #  voxel  # 


Fig.  3.  Numerical  simulations  of  the  SIR-TFE  readout  (a)  and  resulting  qMT  parameter 
fits  for  the  ID  test  object  defined  in  (b-d).  (a)  The  Mzf  for  each  region  (WM  shown 
here)  and  tj  was  simulated  and  reordered  into  the  corresponding  k-space  filter.  The 
resulting  filters  were  applied  to  the  ID  test  object  as  described  in  the  text,  (b-d) 
From  the  simulated  fit  parameters  (solid  gray  lines),  it  can  be  seen  that  the  TFE  readout 
blurs  the  parameter  maps  with  little  or  no  bias  (except  for  kmf  in  CSF  regions,  which  do 
not  exhibit  an  MT  effect). 


the  final  blurred  1 D  object  at  each  tj.  Finally,  to  assess  the  effect  of 
the  TFE  readout  on  qMT  parameter  maps,  the  magnitude  of  the 
blurred  test  object  signal  at  each  voxel  was  fit  to  the  Mzf  component 
of  Eq.  (7)  as  described  in  the  Data  analysis  section. 

Subjects 

MRI  was  performed  on  thirteen  healthy  volunteers  (22-37  years 
old,  10  male,  3  female).  To  test  reproducibility,  six  of  the  healthy 
volunteers  were  asked  to  undergo  a  second  MRI  scan  at  least  two 
weeks  after  the  first  session.  The  study  was  approved  by  our  local  in¬ 
stitutional  review  board,  and  signed  consent  was  obtained  prior  to 
all  examinations. 

Data  acquisition 

Imaging  was  performed  using  a  7.0-T,  Philips  Achieva  MR  scanner 
(Philips  Healthcare,  Best,  The  Netherlands).  A  quadrature  volume  coil 
was  used  for  excitation  and  a  32-channel  head  coil  (Nova  Medical, 
Wilmington,  MA,  USA)  was  used  for  signal  reception.  For  qMT  imag¬ 
ing,  SIR-TFE  data  were  collected  in  each  subject  using  the  general 
pulse  sequence  shown  in  Fig.  1. 

An  initial  experiment  was  performed  in  one  healthy  volunteer  to  de¬ 
termine  the  effect  of  the  post-TFE  saturation  train  and  predelay  time  td 
on  the  qMT  parameter  maps.  For  this  initial  experiment,  SIR-TFE  data 
were  acquired  in  a  single  5-mm  axial  slice  with  and  without  the 
post-TFE  saturation  train  (see  Fig.  1,  shaded  area  labeled  SAT)  over  a 
range  of  td  values  (0.125-10  s).  Additional  imaging  parameters  includ¬ 
ed:  tj  logarithmically  spaced  between  6  ms  and  2  s  (15  values)  and  tj  = 
10  s,  TFE  echoes  per  shot  =  53,  TFE  pulse-to-pulse  interval  (t)/TE/cc  = 
2.8  ms/1.4  ms/15°,  SENSE  factor  =  2,  field-of-view  =  212  x  212  mm2, 
resolution  =  2.0  x  2.0  mm2,  and  number  of  signal  acquisitions  averaged 
(NSA)  =  2. 

Based  upon  the  results  of  this  experiment  along  with  previous  nu¬ 
merical  simulations  (Gochberg  and  Gore,  2007),  a  td  of  2.5  s  was  chosen 
to  balance  the  scan  time  and  SNR  constraints  for  whole-brain  SIR-TFE  im¬ 
aging.  Whole-brain  qMT  data  were  acquired  in  12  volunteers  (six 
scanned  twice)  via  a  three-dimensional  (3D)  SIR-TFE  sequence  using 
the  previously  listed  parameters  except:  tj  logarithmically  spaced  be¬ 
tween  6  ms  and  2  s(13  values)  and  tj  =  8  s,  SENSE  factor  =  4  (2  anteri¬ 
or-posterior,  2  superior-inferior),  field-of-view  =  21 2*  212*  90  mm3, 

resolution  =  2.0 x  2.0 x  3.0  mm3,  and  NSA=1.  This  resulted  in  an  acqui¬ 
sition  time^  19  min  for  30  slices. 

Recall  that  the  signal  model  [Eq.  (7)]  has  terms  (Sfand  Sm)  that  ac¬ 
count  for  the  effect  of  the  inversion  pulse  on  the  free  and  macromo¬ 
lecular  pool  magnetizations.  Sf  was  included  as  a  free  parameter  in 
the  fit  as  described  in  the  Data  analysis  section,  while  Sm  was  numer¬ 
ically  estimated  as  described  in  the  Pulse  sequence  section.  Because 
Sm  is  sensitive  to  Bf  (see  Fig.  2d),  this  numerical  estimation  required 
an  independent  measurement  of  B^“.  Asa  result,  was  estimated  in 
same  volume  as  the  SIR-TFE  data  using  the  actual  flip  angle  imaging 
(AFI)  method  (Yarnykh,  2007)  with  TR1/TR2=  125/25  ms  and  a  60° 
slab-selective  excitation  pulse  (asymmetric  sine  pulse  with  Gaussian 
apodization). 

Data  analysis 

All  data  analyses  were  performed  in  MATLAB  (Mathworks, 
Natick,  MA).  Prior  to  data  fitting,  each  SIR-TFE  and  AFI  volume  was 
co-registered  to  the  SIR-TFE  volume  acquired  at  t,  =  110  ms  (middle 
value)  using  a  3D  rigid  body  registration  based  upon  normalized  mu¬ 
tual  information  (Viola  and  Wells,  1997).  Following  co-registration, 
automatic  brain  extraction  was  performed  (Smith,  2002)  and  qMT 
parameter  maps  were  calculated  in  each  volunteer.  The  SIR-TFE  sig¬ 
nal  model  described  in  Eq.  (7)  has  seven  independent  parameters: 
Rim.  Rif.  Sm,  Sf,  M0f,  PSR  =  M0m/M0f,  and  kmf  (kfm  =  kmfPSR).  As  is  the 
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case  with  pulsed  saturation  methods,  the  signal  dependence  on  R1m 
for  SIR  data  is  weak  (Li  et  al.,  2010).  Therefore,  R1m  was  set  equal  to 
R1ffor  fitting  purposes.  The  parameter  Sm  was  numerically  estimated 
for  each  voxel.  This  required  an  independent  estimate  of  the  actual 
flip  angle  in  each  voxel  (cxactuai),  which  was  calculated  from  the  AFI 
data  using  the  following  relationship  (Yarnykh,  2007): 


%  cos’ 
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where  n  =  TR2/TR1,  r  =  S(TR2)/S(TR2),  and  S  is  the  signal  intensity. 
The  resulting  aactuai  map  was  smoothed  with  a  10x10x9  mm3 
moving-average  filter  to  minimize  the  impact  of  imaging  artifacts. 
Following  this  operation,  the  flip  angle  values  were  converted  to 
factual  values  for  the  composite  inversion  pulse  (see  Fig.  2a)  and 
Sm  was  estimated  using  the  procedure  described  in  the  Pulse 
sequence  section  (see  Fig.  2d).  The  remaining  five  parameters  (R1f,  Sf, 
M0f,  kmfi  and  PSR)  were  estimated  for  each  voxel  by  fitting  SIR-TFE 


data  (14  tj  values)  to  the  component  of  Eq.  (7)  in  a  least-squares 
sense  using  the  procedure  described  in  Dortch  et  al.  (2011). 

SIR-TFE  data  had  a  mean  SNR  per  voxel  of  180  +  50  (range  = 
60-320)  within  the  defined  ROIs,  where  SNR  is  defined  as  M0f  divided 
by  the  standard  deviation  (SD)  of  the  residuals  of  the  fit.  Monte  Carlo 
simulations,  similar  to  those  described  by  Li  et  al.  (2010),  were 
performed  to  predict  the  uncertainty  of  the  fit  parameters  at  these 
SNR  levels.  The  t|  and  td  values  listed  above  were  used  for  these  sim¬ 
ulations.  Additional  simulation  parameters  included:  R-im  =  Rif  = 
0.8  s“  \  kmf  =  15  s“  \  PSR  =  15%,  Sm  =  0.7,  and  Sf=  -  0.95.  Over 
an 

SNR  range  of  60-320,  the  SDs  of  the  fit  PSR,  Rif,  and  kmf  values  were 
0.4-2.2%,  0.01-0.03  s-  \  and  0.9-5. 5  s-  \  respectively.  This 
suggests 

that  PSR  and  R1f  can  be  robustly  determined  from  the  in  vivo  brain 
data  collected  herein.  Consistent  with  previous  studies  (Li  et  al., 
2010),  the  uncertainty  in  kmf  is  expected  to  be  much  larger,  especially 
in  lower  SNR  regions. 

Following  this  fitting  procedure,  qMT  parameter  maps  were 
smoothed  with  a  locally-adaptive  Gaussian  filter  (kernel  size  = 
10x10x9  mm3,  full  width  at  half  maximum  =  1/2  kernel  size)  to 
remove  outliers  that  tended  to  occur  at  tissue  boundaries.  To  per¬ 
form  this  operation,  each  filtered  map  was  subtracted  from  the  raw 
parameter  map,  and  outliers  were  defined  as  voxels  whose  value 
was  three  standard  deviations  above  the  mean  difference  across  all 
voxels.  For  these  outliers,  the  value  in  the  raw  parameter  map  was 
replaced  with  the  value  in  the  filtered  map.  This  process  was  iterated 
until  the  number  of  outliers  was  less  than  the  expected  value  (0.3%  of 
the  total  number  of  voxels). 


Statistics 

Mean  qMT  parameters  (PSR,  R1f,  and  kmf)  were  calculated  within  the 
following  regions-of- interest  (ROI):  head  of  the  caudate,  putamen,  thal¬ 
amus,  genu  and  splenium  of  the  corpus  callosum,  internal  capsule,  coro¬ 
na  radiata,  occipital  WM,  and  frontal  WM.  Statistical  comparisons  were 
performed  on  the  mean  ROI  values  to  evaluate  each  parameter's  i)  vari¬ 
ation  across  ROIs  (i.e.,  regional  differences),  ii)  variation  and  reproduc¬ 
ibility  across  time,  and  iii)  variation  across  volunteers.  To  compare 
parameters  across  WM  regions,  a  non-parametric  Wilcoxon  rank-sum 
test  was  performed,  with  a  pb0.05  deeming  a  significant  difference 
between  ROI  values.  To  evaluate  the  test-retest  reproducibility  of  each 
parameter,  a  Bland-Altman  (BA)  analysis  was  performed.  For  the  BA 
analysis,  the  mean  difference  and  the  limits  of  agreement  (LOA  =  mean 
differenced!  1 .96*  SD)  were  tabulated  across  scans  for  all  ROIs.  Addition¬ 
ally,  a  Wilcoxon  signed-rank  test  was  performed  between  the  test  and 
retest  parameter  values  for  each  ROI,  with  a  p  >  0.05  indicating  a 
non-significant  difference  between  scans  at  each  time  point.  To  assess 
the  test-retest  variability  of  each  parameter  within  eachy^OI,  the  coeffi¬ 
cient  of  variation  was  calculated  from:  CV  S=  M 

(  2b  100,  where 

retest  Vi 

S  is  the  SD  of  the  test-retest  difference  across  subjecjtj|  M  is  the  mean 
value  across  all  test-retest  scans  and  subjects,  and  the  ^  term  accounts 
for  the  propagation  of  uncertainty  from  the  difference  operation.  The 
across-cohort  variability  of  each  parameter  within  each  ROI  was  also 
assessed  via:  CVcoh0rt  =  S/M  *  100,  where  S  is  the  SD  across  the  cohort 
and  M  is  mean  value  across  the  cohort.  All  values  are  reported  as  the 
meandiSD  unless  otherwise  stated. 


Results 

The  results  of  the  numerical  simulations  designed  to  assess  the  effect 
of  TFE  readout  on  qMT  parameter  maps  are  shown  in  Fig.  3.  In  Fig.  3a,  the 
evolution  of  M^  during  the  TFE  readout  is  shown  for  WM  as  a  function  of 
tj.  Note  that  this  evolution  is  related  to  the  k-space  filter  of  the  readout.  It 
can  be  seen  that  the  shape  (width  and  rate  of  decay)  of  the  k-space  filter 
changes  as  a  function  of  tj,  which  manifests  as  a  change  in  object  blurring 
as  a  function  of  tj.  The  effect  of  this  on  the  qMT  parameter  maps  is  shown 
in  Figs.  3b-d.  It  can  be  seen  that  the  resulting  qMT  parameter  maps  are 
smoothed  in  the  phase-encoding  direction  with  little  bias  in  the  fit  pa¬ 
rameters.  It  should  be  noted,  however,  that  PSR  values  were  slightly 
underestimated  in  the  WM  region  of  the  ID  test  object.  Additional  sim¬ 
ulations  indicated  that  this  bias  increased  as  the  size  of  the  WM  region 
decreased. 

Fig.  4  displays  PSR  maps  acquired  with  and  without  application  of 
the  post-TFE  saturation  train  (seethe  SAT  region  in  Fig.  1)  as  a  func¬ 
tion  of  td.  For  scans  with  the  saturation  train,  the  fit  PSR  values  at 
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Fig.  4.  (a)  Maps  of  PSR  as  a  function  of  td  without  (top)  and  with  (bottom)  a  post-TFE  saturation  train  and  (b)  corresponding  mean  (±  SD)  slice-wise  PSR  values.  For  scans  with  the 
saturation  train,  all  PSR  values  were  nearly  identical  to  the  values  at  full  recovery  (dashed  line).  Without  the  saturation  train,  PSR  values  were  increasingly  underestimated  with 
decreasing  td. 
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Fig.  5.  Sample  SIR-TFE  images  (a)  and  model  fit  (b)  from  a  slice  at  the  level  of  the  lateral  ventricles  in  a  healthy  control,  (a)  Images  from  six  of  the  14  inversion  times  are  shown.  Note 
the  characteristic  center  brightening  of  the  images  due  to  Bi*~  inhomogeneities,  (b)  Corresponding  SIR  data  from  a  voxel  in  the  genu  of  the  corpus  callosum.  Note  the  agreement 
between  the  SIR  data  (circles)  and  biexponential  model  [solid  black  line,  Eq.  (7)]  and  the  deviation  from  a  monoexponential  model,  which  is  apparent  at  the  shortest  inversion 
times  shown  in  the  zoomed  inset. 


shorter  td  values  were  nearly  identical  to  those  at  full  recovery  (td  = 
10  s).  Without  this  train,  small  deviations  in  PSR  were  observed  at 
td  =  2.5  s;  and  these  were  more  pronounced  at  td  =  1.25  s.  Thus,  the 
post-TFE  saturation  train  allows  for  reduction  of  td  (and  scan  times) 
with  minimal  parameter  bias. 

Representative  3D  SIR-TFE  data  are  shown  in  Fig.  5.  Fig.  5a  shows  a 
single  slice  at  the  level  of  the  lateral  ventricles  acquired  at  six  of  14  tj 
values.  Note  the  characteristic  center  brightening  due  to  B|  inhomo¬ 
geneities.  Fig.  5b  shows  data  from  a  single  voxel  in  the  genu  of  the 
corpus  callosum  and  the  corresponding  model  fit.  Note  the  agreement 
between  the  SIR-TFE  data  and  the  biexponential  model  described  by 
Eq.  (7).  Additionally,  note  the  deviation  from  monoexponential  re¬ 
covery,  which  is  especially  evident  at  the  shortest  inversion  times. 

Based  upon  these  fits,  maps  of  qMT  parameters  were  generated. 
Recall  that  these  maps  were  filtered  to  reduce  the  impact  of  outliers. 
Fig.  6  displays  representative  qMT  parameter  maps  without  filter¬ 
ing,  with  the  previously  described  locally-adaptive  Gaussian  filter, 
and  with  a  global  Gaussian  filter.  The  locally-adaptive  and  global  fil¬ 
ters  both  removed  outliers  in  the  parameter  maps  (see  arrow  in  the 
top  row  and  the  masks  in  the  bottom  row);  however,  the  locally  fil¬ 
tered  maps  were  blurred  to  a  much  smaller  degree.  As  a  result,  we 
employed  the  locally-adaptive  approach  herein.  For  all  parameter 
maps,  14%  of  all  voxels  in  the  post-brain-extraction  volume  were 
smoothed  using  this  approach.  However,  as  seen  in  the  bottom 
row  of  Fig.  6,  a  majority  of  these  voxels  were  located  along  the 
brain  surface  or  within  the  CSF. 

Fig.  7  displays  results  from  four  representative  slices  in  one 
healthy  subject.  The  qMT  parameters  were  uniform  over  most  of  the 
volume  despite  the  presence  of  large  !1B0  and/or  field  inhomoge¬ 
neities  (as  indicated  by  the  heterogeneity  in  the  Sm  maps).  There 
does,  however,  appear  to  be  some  bias  in  the  qMT  parameter  values 
in  midbrain  slices  (black  arrow),  which  typically  (Moore  et  al., 
2010)  exhibit  the  largest  field  inhomogeneities  and  the  lowest  SNR. 
Nevertheless,  these  data  suggest  that  robust  qMT  parameter  mapping 
can  be  achieved  throughout  most  of  the  brain  using  the  3D  SIR-TFE 
protocol  described  herein. 

ROIs  were  defined  in  a  number  of  WM  and  GM  regions  as  shown 
in  Fig.  8.  The  boxplots  in  the  top  row  of  Fig.  9  display  the  mean  ROI 
qMT  parameters  over  the  12  healthy  volunteers.  For  PSR,  the  mean 
value  across  all  WM  ROIs  (17.6  ±  1 .3%)  was  higher  than  the  values 
across  all  GM  ROIs  (10.3  ±  1.6%).  Additionally,  heterogeneity  within 
WM  PSR  values  was  observed,  but  should  be  interpreted  with  caution 
due  to  the  effect  of  multiple  comparisons.  Nevertheless,  differences  be¬ 
tween  the  following  regions  were  detected:  i)  the  genu  of  the  corpus 
callosum  and  occipital  WM  (p  =  0.026),  ii)  the  genu  of  the  corpus 
callosum  and  the  corona  radiata  (p  =  0.026),  iii)  frontal  and  occipital 


WM  (p  =  0.041),  and  iv)  frontal  WM  and  the  corona  radiata  (p  = 
0.041).  Fit  kmf  values  were  higher  in  GM  (24.4±4.4  s-  1)  than  in  WM 
(14.5  +  1.5  s-  1).  Additionally,  note  the  large,  biased  kmf  values  in  and 


Fig.  6.  Representative  qMT  parameter  (PSR,  kmf,  R1f)  maps  with  and  without  filtering. 
Shown  are  (1st  row)  raw  parameter  maps,  (2nd  row)  parameter  maps  filtered  with 
the  locally-adaptive  Gaussian  filter,  (3rd  row)  parameter  maps  filtered  with  the  global 
Gaussian  filter  (with  an  identical  kernel),  and  (4th  row)  masks  of  the  outliers  detected 
using  the  locally-adaptive  filter.  The  arrow  identifies  a  region  with  biased  PSR  values 
that  are  corrected  by  filtering.  Note  that  the  color-scale  in  these  maps  was  chosen  to 
highlight  the  outliers  and  is  different  than  in  Figs.  4  and  8. 
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Fig.  7.  Representative  parameter  maps  from  one  subject  (four  of  30  slices  are  shown).  The  qMT  parameters  (PSR,  kmf,  R1f)  and  the  inversion  efficiency  Sf  were  uniform  over  most  of 
the  volume  despite  the  presence  of  large  field  inhomogeneities.  There  does,  however,  appear  to  be  some  bias  in  the  qMT  parameter  values  in  midbrain  slices  (black  arrows),  which 
typically  exhibit  the  largest  B0  and  inhomogeneities.  This  results  in  a  deviation  of  Sffrom  —  1  in  these  regions. 


around  areas  containing  CSF,  which  is  likely  a  consequence  of  the  weak 
dependence  of  the  signal  on  kmf  when  PSR  ~  0  (see  the  simulated  data 

in  Fig.  3c).  For  R1f,  differences  between  WM  (0.73  ±  0.03  s-  1)  and  GM 
(0.58 ± 0.05  s-  values  were  also  observed.  The  boxplots  in  Fig.  9  give 
an  indication  of  the  variability  of  each  parameter  across  the  healthy  co¬ 
hort.  To  quantify  this,  the  coefficient  of  variation  was  tabulated  for  each 
ROI,  and  the  mean  value  across  all  ROI  is  given  in  Table  1.  From  this, 
it  can  be  seen  that  the  mean  CVcohort  was  b10%  for  all  of  the 
qMT  parameters,  which  is  not  surprising  given  the  small  age 
range  of  the  healthy  cohort  scanned  herein. 

BA  plots  of  the  observed  difference  in  mean  ROI  qMT  parameters 
between  scans  are  shown  in  the  bottom  row  of  Fig.  9;  and  the  results 
from  this  analysis  are  given  numerically  in  Table  1.  The  mean 


difference  for  all  ROIs  across  scans  was  close  to  zero  for  PSR  (0.0%), 
kmf  (1.2  s-1),  and  R1f  (0.01  s-  1),  indicating  a  lack  of  bias  and  reason¬ 
able  reproducibility.  To  further  test  this,  a  Wilcoxon  signed-rank  test 
was  performed  on  the  test-retest  parameter  values  in  each  ROI.  At 
the  p  =  0.05  level,  no  significant  difference  was  observed  between 
test  and  retest  qMT  parameters  in  any  of  the  ROIs  except  for  kmf  in 
the  genu  of  the  corpus  callosum  (p  =  0.031).  The  test-retest  coeffi¬ 
cient  of  variation  (CVretest)  was  also  tabulated  for  each  metric  to  fur¬ 
ther  assess  each  parameter's  variability  across  time.  As  shown  in 
Table  1,  the  relative  CVretest  values  were  consistent  with  the  corre¬ 
sponding  CVcohort  values,  with  kmf  exhibiting  the  highest  variability. 
In  terms  of  absolute  CV  values,  the  test-retest  variability  was  ap¬ 
proximately  20%  lower  than  the  across  cohort  variability. 


Fig.  8.  Representative  PSR  maps  from  a  single  volunteer  with  corresponding  ROIs  (a  =  corona  radiata;  b  =  occipital  WM;  c  =  frontal  WM;  d  =  corpus  callosum,  genu;  e  =  corpus 
callosum,  splenium;  f  =  internal  capsule;  g  =  head  of  caudate;  h  =  thalamus;  i  =  putamen).  White  and  black  dots  represent  WM  and  GM  ROIs,  respectively.  In  practice,  ROIs  were 
defined  bilaterally  and  results  were  averaged  across  hemispheres.  Here  we  show  ROIs  in  one  hemisphere  for  display  purposes. 
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Fig.  9.  (a-c)  Boxplot  of  the  mean  ROI  qMT  parameters  (he  =  head  of  caudate;  put  =  putamen;  thal  =  thalamus;  owm  =  occipital  WM;  see  =  corpus  callosum,  splenium;  ic  =  in¬ 
ternal  capsule;  cr  =  corona  radiata;  gcc  =  corpus  callosum,  genu;  fwm  =  frontal  WM).  On  each  box,  the  central  mark  is  the  median,  the  edges  of  the  box  are  the  25th  and  75th 
percentiles,  and  the  whiskers  extend  to  the  most  extreme  data  points,  (d-f)  Bland-Altman  plots  of  the  difference  in  parameters  for  WM  (black)  and  GM  (gray)  ROIs  across 
scans.  The  solid  line  is  the  mean  difference,  and  the  dashed  lines  are  the  limits  of  agreement  (mean  difference  ±  1.96  SD). 


Discussion 

This  study  demonstrates  the  feasibility  of  performing  whole-brain 
qMT  measurements  in  the  human  brain  in  vivo  at  high  field.  Pulsed 
saturation  and  SSFP-based  approaches  are  difficult  to  implement  at 
high  field  due  to  RF  power  limitations  and/or  magnetic  field  (BT~ 
and  !1  B0)  inhomogeneities.  In  this  study,  we  employed  the  SIR  qMT 
approach,  which  has  been  suggested  to  be  less  sensitive  to  these  is¬ 
sues.  The  biggest  obstacles  to  overcome  were  i)  the  effect  of  B^and 
!1B0  inhomogeneities  on  the  inversion  pulse  and  the  readout  and  ii) 
the  long  scan  times  associated  with  SIR  imaging.  The  former  of 
these  was  mitigated  by  developing  a  novel  B and  !1<B  insensitive  in¬ 
version  composite  pulse  (Fig.  2)  and  employing  a  low-flip  angle  TFE 
readout;  the  latter  was  mitigated  by  the  efficiency  of  the  TFE  readout 
along  with  additional  protocol  optimization  (e.g.,  reducing  the  num¬ 
ber  of  t,  values  to  14,  applying  SENSE  acceleration  in  two  directions). 
Together  this  resulted  in  a  robust  (Fig.  9),  whole-brain  qMT  imaging 
protocol  with  a  scan  time  of  less  than  20  min. 

Previous  qMT  imaging  studies  at  lower  field  strengths  (Dortch  et 
al.,  2011;  Garcia  et  al.,  2010;  Gloor  et  al.,  2008;  Ropele  et  al.,  2003; 
Sled  and  Pike,  2001;  Sled  et  al.,  2004;  Yarnykh  and  Yuan,  2004) 
have  reported  PSR  values  in  the  range  of  11-16%  and  5-9%  for  WM 
and  GM  structures,  respectively  [PSR  =  F  using  the  notation  of  Sled 
and  Pike  (2000,  2001)  and  M0b  using  the  notation  of  Henkelman  et 
al.  (1993)].  The  PSR  values  presented  herein  (WM:  15-20%,  GM:  9- 


13%)  were  approximately  25%  higher.  PSR  should  be  independent  of 
field  strength,  so  these  differences  may  be  related  to  the  SIR-TFE  se¬ 
quence.  As  previously  discussed,  we  modified  the  inversion  pulse 
and  readout  of  our  3.0-T  SIR-FSE  sequence  to  perform  qMT  imaging 
at  7.0  T.  The  effect  of  the  TFE  readout  on  the  fit  qMT  parameters 
was  assessed  via  numerical  simulations  and  was  found  to  result  in  lit¬ 
tle  bias  in  PSR  (Fig.  3).  However,  it  should  be  noted  that  our  previous 
report  at  3.0  T  employed  a  much  longer  TE  (74  ms)  than  was  employed 
herein  (1.4  ms).  Previous  work  (Bjarnason  et  al.,  2005;  Stanisz  et  al., 
1999)  has  demonstrated  that  MT  contrast  is  TE-dependent  in  WM  due 
to  the  microanatomical  compartition  of  water  into  myelin  and  nonmyelin 
water  spaces.  As  a  result,  it  is  reasonable  to  assume  that  PSR  may  also  ex¬ 
hibit  a  TE-dependence.  In  terms  of  the  inversion  pulse,  we  recognize  that 
PSR  is  sensitive  to  the  macromolecular  pool  lineshape  and  T2m  assump¬ 
tions  used  in  the  numerical  estimation  of  Sm.  Similar  to  our  previous  stud¬ 
ies  (Dortch  et  al.,  2011;  Gochberg  and  Gore,  2007),  we  modeled  the 
macromolecular  pool  using  a  Gaussian  lineshape  (T2m=  10  !Js)  because 
the  Super-Lorentzian  exhibits  an  on-resonance  singularity.  Previous 
work  using  a  1-ms  block  inversion  pulse  at  3.0  T  (Dortch  et  al.,  201 1 ) 
found  that  this  was  a  reasonable  approximation;  however,  this  may 
not  be  true  for  the  longer  (5.5  ms),  higher  power  composite  inver¬ 
sion  pulse  employed  herein.  Additional  work  is  needed  to  explore 
the  field-  and  TE-dependence  of  PSR  values  obtained  via  the  SIR  tech¬ 
nique.  Nevertheless,  the  reported  regional  variation  in  PSR  values  was 
consistent  with  previous  qMT  imaging  studies  (Dortch  et  al.,  2011; 


Table  1 

Test-retest  reproducibility  analysis  of  each  qMT  parameter  (PSR,  R1f,  and  kmf).  Shown  are  the  mean  ±  SD  parameter  values  across  all  ROIs  for  the  test  and  retest  scans,  the  resulting 
mean  paired-difference  between  time-points,  the  limits-of-agreement  (LOA),  and  the  mean  ±  SD  test-retest  coefficient  of  variation  (CVretest)  across  all  ROIs.  For  comparison,  the 
corresponding  across-cohort  coefficient  of  variation  (CVcohort)  is  also  given. 


Parameter 

Test  scan  (mean±SD) 

Retest  scan  (mean  +  SD) 

Difference 

LOA 

CVretest  (%) 

CVcohort  (%) 

PSR  (%) 

15.2+3.9 

15.2+3.7 

0.0 

(-  2.2,  2.1) 

4.9 +  1.5 

5.6  +  1.9 

kmf  (S_1) 

16.8+4.7 

18.0+5.5 

1.2 

(-  3.9,  6.3) 

8.2 +2.4 

9.4  +3.6 

R,f  (s-1) 

0.68+0.08 

0.68+0.08 

0.01 

(-  0.03,  0.04) 

1.9 +  1.4 

3.2 +  1.3 

648 
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Garcia  et  al.,  2010;  Sled  et  al.,  2004;  Underhill  et  al.,  2009);  and  additional 
SIR-TFE  studies  in  bovine  serum  albumin  phantoms  at  7.0  T  (data  not 
shown)  found  a  linear  relationship  between  macromolecular  content 
and  PSR.  Thus,  we  postulate  that  the  regional  differences  in  PSR  values 
reported  herein  are  driven  primarily  by  regional  differences  in  myelin 
content,  although  the  absolute  values  may  be  systematically  larger  than 
reported  by  other  techniques. 

Previous  pulsed  saturation  and  SSFP-based  studies  (Garcia  et  al., 
2010;  Gloor  et  al.,  2008;  Ropele  et  al.,  2003;  Sled  and  Pike,  2001; 
Sled  et  al.,  2004;  Yarnykh  and  Yuan,  2004)  report  kmf  values  [kmf  = 
kf/  F  using  the  notation  of  Sled  and  Pike  (2000,  2001);  kmf  =  R  when 
M0f  =  1  using  the  notation  of  Henkelman  et  al.  (1993)]  in  the  range 
of  20-40  s-  1  across  the  brain.  A  previous  SIR  study  (Dortch  et  al., 
201 1 )  at  3.0  T  reports  kmf  values  that  are  approximately  2-fold  slower 
(10-15  s-  1)  with  values  that  are  slower  in  WM  than  GM,  which  is 
consistent  with  the  results  presented  herein.  The  discrepancies  be¬ 
tween  techniques  are  not  surprising  given  the  reported  difficulty  of 
using  pulsed  saturation  to  determine  kmf  (Portnoy  and  Stanisz, 
2007).  In  terms  of  the  current  study,  it  should  be  noted  that  kmf 
showed  the  largest  variability  of  the  qMT  parameters,  which  is  con¬ 
sistent  with  the  results  from  the  Monte  Carlo  simulations.  We  do 
not  expect  this  to  be  a  significant  drawback  as  kmf  has  been  shown 
to  be  insensitive  to  the  pathological  changes  in  spinal  cord  WM  (Smith 
et  al.,  2009). 

While  there  have  been  no  previous  reports  of  R1f  in  human  brain 
at  7.0  T,  it  can  be  shown  that  the  observed  T-,  typically  reported 
\s~  1  /Rif.  Using  this  relationship,  the  mean  WM  and  GM  observed 
Ti  values  were  1372  and  1724  ms,  respectively,  which  are  within 
the  range  of  previously  reported  values  in  human  brain  at  7.0  T 
(Wright  et  al.,  2008).  As  expected,  we  noted  a  significant  correlation 
between  R1f  and  PSR  in  the  healthy  human  brain  (data  not  shown); 
however,  T1  is  sensitive  to  overall  tissue  composition  [e.g.,  water  con¬ 
tent  (Kiricuta  and  Simplaceanu,  1975)]  and  is  believed  to  be  a  less 
specific  marker  for  myelin  in  WM. 

The  increased  SNR  available  at  7.0  T  was  used  here  to  decrease 
scan  time  (~  40  s/slice)  and  increase  resolution  (2x2x3  mm3)  rela¬ 
tive  to  our  3.0-T  protocol.  Moving  forward,  it  may  be  advantageous  to 
look  at  higher  resolution  protocols.  If  we  assume  that  all  imaging  pa¬ 
rameters  are  the  same,  increasing  the  resolution  to  1x1x3  mm3 
would  result  in  an  approximately  two-fold  decrease  in  SNR  (~  70  at 
thermal  equilibrium,  assuming  we  increase  the  number  of  acquired 
points  to  hold  the  field-of-view  constant).  Based  upon  previous  sim¬ 
ulation  work  (Li  et  al.,  2010)  as  well  as  the  simulation  work  presented 
herein,  this  would  be  sufficient  to  robustly  fit  qMT  parameters  over 
most  of  the  brain.  Thus,  it  appears  that  high-resolution  qMT  imaging 
may  be  feasible  in  the  human  brain  in  vivo  at  7.0  T  using  a  protocol 
similar  to  that  described  herein. 

Conclusions 

The  results  of  this  study  demonstrate  the  feasibility  of  performing 
qMT  imaging  in  human  brain  in  vivo  at  high  field.  The  developed 
SIR-TFE  protocol  allowed  for  whole-brain  qMT  imaging  in  less  than 
20  min.  In  healthy  subjects,  intra-subject  reliability  (i.e.,  test-retest) 
was  demonstrated  despite  large  !1B0and  Bf*"  variations.  Additionally,  a 
high  level  of  inter-subject  reproducibility  was  demonstrated  for  the 
qMT  parameters.  Future  work  includes  investigating  high-resolution 
protocols  to  look  at  cortical  features  of  qMT  parameters  and  application 
of  the  approach  in  a  cohort  of  multiple  sclerosis  patients. 
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